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We review the BCS to Bose Einstein condensation (EEC) crossover scenario which is based on the well known 
crossover generalization of the BCS ground state wavefunction ^o- While this ground state has been summarized 
extensively in the literature, this Review is devoted to less widely discussed issues: understanding the effects of 
finite temperature, primarily below Tc, in a manner consistent with ^q. Our emphasis is on the intersection 
of two important problems: high Tc superconductivity and superfluidity in ultracold fermionic atomic gases. 
We address the "pseudogap state" in the copper oxide superconductors from the vantage point of a BCS-BEC 
crossover scenario, although there is no consensus on the applicability of this scheme to high Tc- We argue that 
it also provides a useful basis for studying atomic gases near the unitary scattering regime; they are most likely 
in the counterpart pseudogap phase. That is, superconductivity takes place out of a non-Fermi hquid state where 
preformed, metastable fermion pairs are present at the onset of their Bose condensation. As a microscopic basis 
for this work, we summarize a variety of T-matrix approaches, and assess their theoretical consistency. A close 
connection with conventional superconducting fluctuation theories is emphasized and exploited. 
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I. INTRODUCTION TO QUALITATIVE CROSSOVER 
PICTURE 

This review addresses an extended form of the BCS theory 
of superfluidity (or superconductivity) known as "BCS-Bose 
Einstein condensation (EEC) crossover theory". Here we con- 
template a generalization of the standard mean field descrip- 
tion of superfluids now modified so that one allows attractive 
interactions of arbitrary strength. In this way one describes 
how the system smoothly goes from being a superfluid sys- 
tem of the BCS type (where the attraction is arbitrarily weak) 
to a Rose Einstein condensate of diatomic molecules (where 
the attraction is arbitrarily strong). One, thereby, sees that 
BCS theory is intimately connected to EEC. Attractive inter- 
actions between fermions are needed to form "bosonic-like" 
molecules (called Cooper pairs) which then are driven statisti- 
cally to Eose condense. ECS theory is a special case in which 
the condensation and pair formation temperatures coincide. 

The importance of obtaining a generalization of ECS theory 
which addresses the crossover from ECS to EEC at general 
temperatures T < Tc cannot be overestimated. ECS theory 
as originally postulated can be viewed as a paradigm among 
theories of condensed matter systems; it is complete, in many 
ways generic and model independent, and well verified ex- 
perimentally. The observation that a ECS-like approach goes 
beyond strict ECS theory, suggests that there is a larger mean 
field theory to be addressed. Equally exciting is the possibil- 
ity that this mean field theory can be discovered and simul- 
taneously tested in a very controlled fashion using ultracold 
fermionic atoms. We presume here that it may also have ap- 
plicability to other short coherence length materials, such as 
the high temperature superconductors. 

This Review is written in order to convey important new 
developments from one subfield of physics to another. Here it 
is hoped that we will communicate the recent excitement felt 
by the cold atom community to condensed matter physicists. 
Conversely we wish to communicate a comparable excitement 
in studies of high temperature superconductivity to atomic 
physicists. There has been a revolution in our understanding 
of ultracold fermions in traps in the last few years. The mile- 
stones in this research were the creation of a degenerate Fermi 
gas (1999), the formation of dimers of fermions (2003), Eose- 
Einstein condensation (EEC) of these dimers (late 2003) and 
finally superfluidity of fermionic pairs (2004). 

In high temperature superconductors the ECS-EEC 
crossover picture has been investigated for many years now. It 
leads to a particular interpretation of a fascinating, but not well 
understood phase, known as the "pseudogap state". In the cold 
atom system, this crossover description is not just a scenario, 
but has been realized in the laboratory. This is because one 
has the ability (via Feshbach resonances) to apply magnetic 
fields in a controlled way to tune the strength of the attrac- 
tion. We note, finally, that research in this field has not been 
limited exclusively to the two communities. One has seen 
the application of these crossover ideas to studies of excitons 
in so lids (Lai et ai, 2004), to nuclear physics (Eald^^^^, 
[1995; Heiselbere, 2004a) and to part icle physics ([i takun , 
l2003^.KIinkhamer and Volovik..2004:.Kolb and Heinzll2003 ) 



as well. There are not many problems in physics which have 
as great an overlap with different subfield communities as this 
"ECS-EEC crossover problem". 

In this Review, we begin with an introduction to the qual- 
itative picture of the ECS-EEC crossover scenario which is 
represented schematically in Fig. [2 We discuss this picture 
in the context of both high Tc superconductivity and ultra- 
cold atomic superfluidity, summarizing experiments in both. 
In Section II, we present a more quantitative study of this 
crossover primarily at and below the superfluid transition tem- 
perature Tc- Section III contains a detailed discussion of the 
superfluid density and gauge invariance issues, in large part 
to emphasize self consistency checks for different theoreti- 
cal approaches to crossover. These alternative theories are 
summarized and compared in Section IV. Sections V and VI, 
respectively, present the physical implications of crossover 
physics, in one particular rendition, in the context of ultracold 
fermionic superfluids and high temperature superconductors. 
Our conclusions are presented in Section VII and a number 
of Appendices are included for clarity. A reader interested in 
only the broad overview may choose to omit Sections II-IV. 

The convention and notation as well as a list of abbrevia- 
tions used in this Review are listed in Appendix IDlat the end. 



A. Fermionic Pseudogaps and Meta-stable Pairs: Two 
Sides of the Same Coin 

A number of years ago Eagles llEaglesLll969i) and Leggett 
( iLegge tt. 1980) independently noted' that the ECS ground 
state wavefunction 

*o = nkK + i^k4,icLfc,i)|o) (1) 

had a greater applicability than had been appreciated at the 
time of its original proposal by Eardeen, Cooper and Schrief- 
fer (ECS). As the strength of the attractive pairing interaction 
U {< 0) between fermions is increased, this wavefunction is 
also capable of describing a continuous evolution from ECS 
like behavior to a form of Eose Einstein condensation. What 
is essential is that the chemical potential /i of the fermions be 
self consistently computed as U varies. 

The variational parameters Uk and iik which appear in Eq.^ 
are usually represented by the two more directly accessible 
parameters Asc(O) and /i, which characterize the fermionic 
system. Here Asc(O) is the zero temperature superconduct- 
ing order parameter These fermionic parameters are uniquely 
determined in terms of U and the fermionic density n. The 
variationally determined self consistency conditions are given 
by two ECS-like equations which we refer to as the "gap" and 



Through private correspondence it appears that F. Dyson may have made 
similar observations even earlier. 
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BCS 

weak coupling 

large pair size 
k- space pairing 

strongly overlapping 
Cooper pairs 



BEC 

strong coupling 

small pair size 
r-space pairing 

ideal gas of 
preformed pairs 




T* = T, 



T* > 



FIG. 1 Contrast between BCS- and BEC-based superfluids. Here 
T* represents the temperature at which pairs form, while Tc is that 
at which they condense. 



"number" equations respectively. 



2Ek 



1 - 



where ^ 



(2) 
(3) 

(4) 



and = k"^ /2m is the fermion energy dispersion. Within 
this ground state there have been extensive studies of collec- 
tive modes (Cote and Griffin, 1993; Micnas et ai, 1990; Ran- 
deria. . 19951 effects of two dimensionality ( Randeria,^9950, 
and, more recently, extensions to atomic gases (Peral ief al\. 
|2003; Viverit et al., 2004). Nozieres and Schmitt-Rink were 
the first (Nozier e s and Schmitt-Rink. J_98 5) to address nonzero 
T. We will outline some of their conclusions later in this Re- 
view. Randeria and co-workers reformulated the approach of 
Nozieres and Schmitt-Rink (NSR) and moreover, raised the 
interesting possibility that crossover ph ysics might be r elevant 
to high temperature superconductors ( Randerialll995l) . Sub- 
sequently other workers have applied this picture to the high 



^ In general, the order parameter Asc and the excitation gap may have a 
k-dependence, Aj. = Ai/Ji^, except for a contact potential. To keep 
the mathematical expressions simpler, however, we neglect i/p^ for the 
s-wave pairing under investigation of this Review. It can be easily rein- 
serted hack when needed in actual calculations. In particular, we use 
fk = exp(— fc^/2fc^) with a large cutoff fcc in most of our numerical 
calculations. 

^ Throughout this Review, we set 7l = 1 and kg = 1. 
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FIG. 2 Typical behavior of the T — Q chemical potential in the three 
regimes. As the pairing strength increases from 0, the chemical po- 
tential ^ starts to decrease and then becomes negative. The character 
of the system changes from fermionic (/i > 0) to bosonic (/i < 0). 
The PG (pseudogap) case corresponds to non-Fermi liquid based su- 
perconductivity, and Uc corresponds to critical coupling for forming 
a two fermion bound state in vacuum, as in Eq. jl7> below. 



Tc cuprates ("Chen et al}. 


'1998*; "Micnas and Robaszkiewiczl 


.1998: Ranninser and Robin. 


1996) and ultracold fermions 


jMilstein et adl2002; Ohashi and Griffin. 2002) as well as for- 



mulated alternative schem es (.Griffin and Ohashi, 2003; Fieri 
et a/.^004) for addressing T ^ 0. Importantly, a number of 
experimentalists, most notably Uemura (Uemura, 1997 ), have 



.1999: 



claimed evi dence in support ( Deutsche^ 

11999; Renn er ef ad 11998) of the BCS-BEC crossover picture 
for high Tc materials. 

Compared to work on the ground state, considerably less 
has been written on crossover effects at nonzero temperature 
based on Eq. (^. Because our understanding has increased 
substantially since the pioneering work of NSR, and because 
they are the most interesting, this review is focussed on these 
finite T effects, at the level of a simple mean field theory. 
Mean field approaches are always approximate. We can as- 
cribe the simplicity and precision of BCS theory to the fact 
that in conventional superconductors the coherence length ^ is 
extremely long. As a result, the kind of averaging procedure 
implicit in mean field theory becomes nearly exact. Once ^ 
becomes small, BCS is not expected to work at the same level 
of precision. Nevertheless even when they are not exact, mean 
field approaches are excellent ways of building up intuition. 
And further progress is not likely to be made without investi- 
gating first the simplest of mean field approaches, associated 
with Eq. (0. 

The effects of BEC-BCS crossover are most directly re- 
flected in the behavior of the fermionic chemical potential fi. 
We plot the behavior of ji in Fig.|2j which indicates the BCS 
and BEC regimes. In the weak coupling regime fi ~ Ep and 
ordinary BCS theory results. However at sufficiently strong 
coupling, /i begins to decrease, eventually crossing zero and 
then ultimately becoming negative in the BEC regime, with 
increasing \U\. We generally view = as a crossing point. 
For positive fi the system has a remnant Fermi surface, and we 
say that it is "fermionic". For negative fi, the Fermi surface is 
gone and the material is "bosonic". 

The new and largely unexplored physics of this problem 
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lies in the fact that once outside the BCS regime, but before 
BEC, superconductivity or superfluidity emerge out of a very 
exotic, non-Fermi liquid normal state. Indeed, it should be 
clear that Fermi liquid theory must break down somewhere 
in the intermediate regime, as one goes continuously from 
fermionic to bosonic statistics. This intermediate case has 
many names: "resonant superfluidity", the "unitary regime" 
or, as we prefer (because it is more descriptive of the many 
body physics) the "pseudogap (PG) regime". The PG regime 
plotted in Fig.|2lis drawn so that its boundary on the left co- 
incides with the onset of an appreciable excitation gap at Tc, 
called A(Tc); its boundary on the right coincides with where 
/J, = 0. This non-Ferm i liquid bas ed superconducti vity has 
been extensively studied llChen et alV [1298, 2000; Jank o et a/.L 
ll997l:lMalv^rfl?lll999ah ^ 

Fermi liquid theory is inappropriate here principally be- 
cause there is a gap for creating fermionic excitations. The 
non-Fermi liquid like behavior is reflected in that of the 
fermionic spectral function (.Malv et ai. 1999a) which has a 
two-peak structure in the normal state, broadened somewhat, 
but rather similar to its counterpart in BCS theory. Impor- 
tantly, the onset of superconductivity occurs in the presence of 
fermion pairs, which reflect this excitation gap. Unlike their 
counterparts in the BEC limit, these pairs are not infinitely 
long lived. Their presence is apparent even in the normal 
state where an energy must be applied to create fermionic ex- 
citations. This energy cost derives from the breaking of the 
metastable pairs. Thus we say that there is a "pseudogap" at 
and above Tc- But pseudogap effects are not restricted to this 
normal state. Pre-formed pairs have a natural counterpart in 
the superfluid phase, as noncondensed pairs. 

It will be stressed throughout this Review that gaps in the 
fermionic spectrum and bosonic degrees of freedom are two 
sides of the same coin. A particularly important observation 
to make is that the starting point for crossover physics is based 
on the fermionic degrees of freedom. Bosonic degrees of free- 
dom are deduced from these; they are not primary. A nonzero 
value of the excitation gap A is equivalent to the presence of 
metastable or stable fermion pairs. And it is only in this indi- 
rect fashion that we can probe the presence of these "bosons", 
within the framework of Eq. Q. 

In many ways this crossover theory appears to represent a 
more generic form of superfluidity. Without doing any calcu- 
lations we can anticipate some of the effects of finite tempera- 
ture. Except for very weak coupling, pairs form and condense 
at different temperatures. The BCS limit might be viewed as 
the anomaly. Because the attractive interaction is presumed 
to be arbitrarily weak, in BCS the normal state is unaffected 
by U and superfluidity appears precipitously, that is without 
warning at Tc. More generally, in the presence of a moder- 
ately strong attractive interaction it pays energetically to take 
some advantage and to form pairs (say roughly at temperature 
T*) within the normal state. Then, for statistical reasons these 
bosonic degrees of freedom ultimately are driven to condense 
atTc < T*,just as in BEC. 

Just as there is a distinction between Tc and T*, there must 
be a distinction between the superconducting order parame- 
ter A^c and the excitation gap A. The presence of a normal 




\ T* 

FIG. 3 Contrasting behavior of the excitation gap A(r) and super- 
fluid order parameter AsciT) versus temperature. The height of the 
shaded region roughly reflects the density of noncondensed pairs at 
each temperature. 



state excitation gap or pseudogap for fermions is inextricably 
connected to this generalized BCS wavefunction. In Figure|3] 
we present a schematic plot of these two energy parameters. 
It may be seen that the order parameter vanishes at Tc, as in 
a second order phase transition, while the excitation gap turns 
on smoothly below T* . It should also be stressed that there is 
only one gap energy scale in the ground state ( Leggett. 198Ci) 
of Eq. Q. Thus at zero temperature 

A,e(0) = A(G) . (5) 

In addition to the distinction between A and A^c, another 
important way in which bosonic degrees of freedom are re- 
vealed is indirectly through the temperature dependence of A. 
In the BEC regime where fermion pairs are pre-formed, A is 
essentially constant for all T < Tc (as is fj,). By contrast in 
the BCS regime it exhibits the well known temperature depen- 
dence of the superconducting order parameter. This is equiva- 
lent to the statement that bosonic degrees of freedom are only 
present in the condensate for this latter case. This behavior is 
illustrated in Fig.|3 

Again, without doing any calculations we can make one 
more inference about the nature of crossover physics at finite 
T. The excitations of the system must smoothly evolve from 
fermionic in the BCS regime to bosonic in the BEC regime. 
In the intermediate case, the excitations are a mix of fermions 
and meta-stable pairs. Figure |5] characterizes the excitations 
out of the condensate as well as in the normal phase. This 
schematic figure will play an important role in our thinking 
throughout this review. In the BCS and BEC regimes one is 
led to an excitation spectrum with a single component. For the 
PG case it is clear that the excitations now have two distinct 
(albeit, strongly interacting) components. 



B. Introduction to high Tc Superconductivity: Pseudogap 
Effects 

This Review deals with the intersection of two fields and 
two important problems: high temperature superconductors 
and ultracold fermionic atoms in which, through Feshbach 
resonance effects, the attractive interaction may be arbitrarily 
tuned by a magnetic field. Our focus is on the broken symme- 
try phase and how it evolves from the well known ground state 
at r = to T = Tc- We begin with a brief overview (Tallon 
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T 



FIG. 4 Comparison of typical temperature dependences of the exci- 
tation gaps in tlie BCS (left) and BEC (right) limits. For the former, 
the gap is small and vanishes at Tc\ whereas for the latter, the gap is 
very large and essentially temperature independent. 



BCS 



Pseudogap (PG) 



BEC 




FIG. 5 



The 



character of the excitations in the BCS-BEC crossover 



both abwe md below Tc. The excitations are primarily fermionic 
BogoliuTiov quasiparticles in the BCS limit and bosonic pairs (or 
"Feshbcch bosons") in the BEC limit. In the PG case the "virtual 
molecu es" ■ ;onsist primarily of "Cooper" pairs of fermionic atoms. 



and Loram, l200lL iTimusk and Statd 1 19991) of pseudogap ef- 
fects in the hole doped high temperature superconductors. A 
study of concrete data in these systems provides a rather nat- 
ural way of building intuition about non-Fermi liquid based 
superfluidity, and this should, in turn, be useful for the cold 
atom community. 

It has been argued by some ( Chen et al. 1 999 ; Micnas and 
Robaszkiewi cz,.1998; Pieri and Strinati. 200 0; Ranninger and 
Robin. ll996l:lYanase et adl2003i) fliat aBCS-BEC crossover- 
induced pseudogap is the origin of the mysterious normal state 
gap observed in high temperature superconductors. Although, 
this is a highly contentious subject, a principal supporting ra- 
tionale is that this above-Tc excitation gap seems to evolve 
smoothly into the fermionic gap within the superconducting 
state. This, in conjunction with a panoply of normal state 
anomalies which appear to represent precursor superconduc- 
tivity effects, suggests that the pseudogap has more in com- 
mon with the superconducting phase than with the insulating 
magnetic phase associated with the parent compound. 

The arguments in favor of this (albeit, controversial) 
crossover viewpoint rest on the following observations: (i) the 
transition temperature is anomalously high and the coherence 
length ^ for superconductivity anomalously short, around 
lOA as compared more typically with lOOOA. (ii) These sys- 
tems are quasi-two dimensional so that one might expect pre- 
formed pairs or precursor superconductivity effects to be en- 
hanced, (iii) The normal state gap has the same d-wave sym- 
metry as the superconducting order parameter (Campuzano 
et al, 2004; Damascelli et al, 2003), and (iv) to a good ap- 
proximation its onset temperature, called T*, satisfies (Ku- 



ns^^VCm'/rrie) or -^Gn^iJ i^' / "^e) (cm ^) 
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FIG. 6 Uem ura plot iUemural 



tion temperature as a function t 



llll 



1997^ of the superconducting transi- 



f Fermi temperature (lower axis) and 



of the low T superfluid density (upper axis) for various supercon- 



ductors. Tme plot indicates ho\ ' "exotic" superconductors appear to 
be in a distinct category as compared with more conventional metal 
supercondijCtors such as Nb, 5 n, Al and Zn. Note the logarithmic 
scales. 



gler et al., l200lHOda eT7l. 



. 11997b T* w 2A(0)/4.3 which 
satisfies the BCS scaling relation, (v) There is widespread 
evidence for pseudogap effects both_aboye_(T allon and Lo- 
ram, 2001; Timusk and Statt, 1999) and below ( iLoram et al\. 
1994; Stajic et al, 2003b) Tc. Finally, it has also been argued 
that short coherence length superconductors may qu ite gener- 
ally e xhibit a distinctive form of superconductivity (lUemural 
Il997b . Figure |5] shows a plot of data collected by Uemura 
which seems to suggest that the traditional BCS superconduc- 
tors stand apart from other more exotic (and frequently short 
forms of superconductivity. From this plot one can infer, 
that, except for high Tc, there is nothing special about the high 
Tc superconductors; they are not alone in this distinctive class 
which includes the organics and heavy fermion superconduc- 
tors as well. Thus, to understand them, one might want to 
focus on this simplest aspect (short ^) of high Tc superconduc- 
tors, rather than invoke more exotic and less generic features. 

In Fig. Q we show a sketch of the phase diagram for the 
copper oxide superconductors. Here x represents the concen- 
tration of holes which can be controlled by adding Sr substi- 
tutionally, say, to Lai_2;Sr2;Cu04. At zero and small x the 
system is an antifeiTomagnetic (AFM) insulator. Precisely at 
half filling {x = 0) we understand this insulator to derive from 
Mott effects. These Mott effects may or may not be the source 
of the other exotic phases indicated in the diagram, as "SC" 
for superconductivity and the pseudogap phase. Once AFM 
order disappears the system remains insulating until a criti- 
cal concentration (typically around a few percent Sr) when 
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Hole concencentration x 

FIG. 7 Typical phase diagram of hole-doped high Tc superconduc- 
tors. There exists a pseudogap phase above Tc in the underdoped 
regime. Here SC denotes superconductor, and T* is the temperature 
at which the pseudogap smoothly turns on. 



A(T^ ) 600 




0.00 0.05 0.10 0.15 0.20 

hole concentration x 

FIG. 8 Pseudogap magnitude at Tc versus doping concentration in 
hole doped high Tc supe rconductors measured wi th different experi- 
mental techniques ffrom lTallon and Loramll2boih . Units on vertical 
axis are temperature in degrees K. 




150 200 

T(K) 

FIG. 9 (Color) Temperature dependence of the excitation gap from 
ARPES measurements for optimally doped (filled black circles), un- 
derdoped (red squares) and highly underdoped (blue inverted trian- 
gles) single-crystal BSCCO samples (taken from Dine et al, 199g). 
There exists a pseudogap phase above Tc in the underdoped regime. 
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Sample 



100 
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300 



FIG. 10 STM measurements of the dl /dV characteristics of an un- 
derdoped BSCCO sample inside (solid) and outside (dashed) a vortex 

Il ll998l dl/dV is proportional to the 



core at low T from lRenner et 
density of states. 



an insulator-superconductor transition is encountered. Here 
photoemission studies (Campuzano et al. 2004; Damascelli 
et al, 2003) suggest that once this line is crossed, /i appears 
to be positive. For x < 0.2, the superconducting phase has a 
non-FOTnjJigui^^gr^seu^ggaDDed) normal state (Tallon and 
Loram7 l200ll) . We note an important aspect of this phase di- 
agram at low X. As the pseudogap becomes stronger (T* in- 
creases), superconductivity as reflected in the magnitude of 
Tc becomes weaker. One way to think about this competi- 
tion is through the effects of the pseudogap on Tc- As this 
gap increases, the density of fermionic states at Ep decreases, 
so that there are fewer fermions around to participate in the 
superconductivity. 

Importantly, there is a clear excitation gap present at the 
onset of superconductivity which is called A(Tc) and which 
vanishes at x ^ 0.2. Th e magnitude of the pseud ogap at Tc 
is shown in Fig.|H| from (Tal lon and Lorarrll200ll) . Quite re- 
markably, as indicated in the figure, a host of different probes 
seem to converge on the size of this gap. 

Figure|9]indicates the temperature dependence of the exci- 
tation gap for three different hole stoichiometrics. These data 



were taken llDingef a/.LIl996h from angle resolved photoemis- 
sion spectroscopy (ARPES) measurement on single-crystal 
Bi2Sr2CaCu208+5 (BSCCO) samples. For one sample shown 
as circles, (coiTesponding roughly to "optimal" doping) the 
gap vanishes roughly at Tc as might be expected for a BCS 
superconductor At the other extreme are the data indicated 
by inverted triangles in which an excitation gap appears to be 
present up to room temperature, with very little temperature 
dependence. This is what is referred to as a highly underdoped 
sample (small x), which from the phase diagram can be seen 
to have a rather low Tc- Moreover, Tc is not evident in these 
data on underdoped samples. This is a very remarkable fea- 
ture which indicates that from the fermionic perspective there 
appears to be no profound sensitivity to the onset of supercon- 
ductivity. 

While the high Tc community has focussed on pseudogap 
effects above Tc, there is a good case to be made that these ef- 
fects also persist beloWj_ShowninFi 2 . 1 1 Olare STM data (Ren- 
ner et a/.7 ri998,) taken below Tc within a vortex core (solid 
lines) and between vortices in the bulk (dashed lines). The 
quantity dl/dV may be viewed as a measure of the fermionic 
density of states at energy E given by the voltage V. This 
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figure shows that there is a clear depletion of the density of 
states around the Fermi energy {V = 0) in the normal phase 
within the core. Indeed the size of the inferred energy gap 
(or pseudogap) corresponds to the energy separation between 
the two maxima in dl/dV; this can be seen to be the same 
for both the normal and superconducting regions (solid and 
dashed curves). This figure emphasizes the fact that the ex- 
istence of an energy gap has little or nothing to do with the 
existence of phase coherent superconductivity. It also sug- 
gests that pseudogap effects effectively persist below T^, the 
normal phase underlying superconductivity for T < Tc is not 

a Fermi liquid. 

Analysis of thermodvnamical data dLoram et a/.LIl994t Tal- 
lon and Loram, 2001) has led to a similar inference. Figure fTTI 
presents a schematic plot of the entropy S and specific heat 
for the case of a BCS superconductor, as contrasted with a 
pseudogapped superconductor. Actual data are presented in 
Fig- ED Figure ^2 makes it clear that in a BCS supercon- 
ductor, the normal state which underlies the superconducting 
phase, is a Fermi liquid; the entropy at high temperatures ex- 
trapolates into a physically meaningful entropy below Tc- For 
the PG case, the Fermi liquid-extrapolated ent ropy becomes 
negat ive. In this way Loram and co-workers dLoram et ali 
11994 deduced that the normal phase underlying the super- 
conducting state is not a Fermi liquid. Rather, they claimed to 
obtain proper thermodynamics, it must be assumed that this 
state contains a persistent pseudogap. They argued for a dis- 
tinction between the excitation gap A and the superconduct- 
ing order parameter, within the superconducting phase. To fit 
their data they presumed a modified fermionic dispersion 



where 



(6) 



(7) 



Here Apg is taken on phenomenological grounds to be T- 
independent. The authors argue that the pseudogap contribu- 



tion may arise fr om physics unrelate d to the superconductiv- 
ity. While others fchakravarty et al. , 2001; Nozieres and Pis- 
tolesi, 1999) have similarly postulated that Apg may in fact 
derive from another (hidden) order parameter or underlying 
band structure effect, in general the fermionic dispersion rela- 
tion £'k will take on a different character from that assumed 
above, which is very specific to a superconducting origin for 
the pseudogap. 

Finally, Fig. [2] makes the claim for a persistent pseudo- 
gap below Tc in an even more suggestive way. Figure fT^ a) 
represents a schematic plot of excitation gap data such as are 
shown in Fig.|9] Here the focus is on temperatures below Tc- 
Most importantly, this figure indicates that the T dependence 
in A varies dramatically as the stoichiometry changes. Thus, 
in the extreme underdoped regime, where PG effects are most 
intense, there is very little T dependence in A below Tc- By 
contrast at high x, when PG effects are less important, the 
behavior of A follows that of BCS theory. What is most im- 
pressive however, is that these wide variations in A(T) are 
not reflected in the superfluid density ps (T), which is propor- 
tional to the inverse square of the London magnetic penetra- 
tion depth and can be measured by, for example, muon spin 
relaxation (/iSR) experim ents. Necessarily, Ps{T) vanishes 
at Tc. What is plotted (iStaiic et al\. l2003bh in Fig. Ob) is 
Ps{T) — Ps{0) versus T. That these data all seem to sit on a 
rather universal curve is a key point. The envelope curve in 
Ps{T) — Ps(0) is associated with an "optimal" sample where 
A(T) essentially follows the BCS prediction. Figure [T2l then 
indicates that, despite the highly non-universal behavior for 
A(T), the superfluid density does not make large excursions 
from its BCS- predicted form. This is difficult to understand if 
the fermionic degrees of freedom through A(T) are dominat- 
ing at all X. Rather this figure suggests that something other 
than fermionic excitations is responsible for the disappearance 
of superconductivity, particularly in the regime where A(T) 
is relatively constant in T. At the very least pseudogap effects 
must persist below Tc. 

Driving the superconductivity away is another important 
way to probe the pseudogap. This occurs naturally with 
temperatures in excess of Tc, but it also occur s when suffi- 
cient pair breaking is present throug h impur i ties ([Tallon et alX. 
ll997l: lTallon and Loram, 2001; Zha o ef adri993h or applied 
magnetic fields dAndo et ali 119951) . An important effect of 
temperature needs to be stressed. With increasing T > Tc, 
the d-wave shape of the excitation gap is rapidly washed out 
dCampu zano et aZ.V2004': ^amascelli et a/.L 12003'); the nodes 
of the order parameter, in effect, begin to spread out, just 
above Tc. The inverse of this effect should also be empha- 
sized: when approached from above, Tc is marked by the 
abrupt onset of long lived quasiparticles. 

One frequently, and possibly universally, sees a 
superconductor-insulator (SI) transition when Tc is driven 
to zero in the presence of a pseudogap. This suggests a 
simple scenario: that the pseudogap may survive when 
superconductivity is suppressed. In this way the ground state 
is no longer a simple metal. This SI transition is see n upon 
Zn doping ( Tallon et al. , 1997; Tallon and Loram, 2001 ; Zhao 
et al. .1993 ). as well as in the presence of applied magnetic 
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FIG. 12 Temperature dependence of ( a) the fermionic excitation gap A and (b) superfluid density ps from experiment for various doping 
concentrations (from iStaiic ef~ Figure (a) siiows a schematic plot of A as seen in the ARPES and STM tunneling data. When 

A(rc) 7^ 0, there is little correlation between A(T') and Ps{T); this figure suggests that something other than fermionic quasiparticles (e.g., 
bosonic excitations) may be responsible for the disappearance of superconductivity with increasing T. Figure (b) shows a quasi-universal 
behavior for the slope dps/dT at different doping concentrations, despite the highly non-universal behavior for A(T). 



fields (lAndo et a/.LIl995l) . Moreover, the intrinsic change in 
stoichiometry illustrated in the phase diagram of Fig. also 
leads to an SI transition. One can thus deduce that the effects 
of pair-breaking on Tc and T* are very different, with the 
former being far more sensitive than the latter. 

It is clear that pseudogapped fermions are in evidence in a 
multiplicity of experiments. Are there indications of bosonic 
degrees of freedom in the normal state of high Tc super- 
conductors? The answer is unequivocally yes: meta-stable 
bosons are observable as superconducting fluctuations. These 
bosons are intimately connected to the pseudogap; in effect 
they are the other side of the same coin. These bosonic ef- 
fects are enhanced in the presence of the quasi-two dimen- 
sio nal lattice structure of these m aterials. Very detailed analy- 
ses ( iLarkin and Varlamovll200ll) of thermodynamic and trans- 
port properties of the highest Tc or "optimal" samples reveal 
clea rly these pre-formed pair s. Moreover they are responsi- 
ble jKwon and Dorsevlll999h for divergences at Tc in the dc 
conductivitv a and ii^i^^nsverse thermoelectric (Ussishkin 
et ai, 2002) response axy. These transport coefficients are 
defined more generally in terms of the electrical and heat cur- 
rents by 



jelec ^ crE + a( - VT) 



rheat 



(iE + k( - VT) 



(9) 



Here a is the conductivity tensor, k the thermal conductivity 
tensor, and a, a are thermoelectric tensors. Other coefficients, 
while not divergent, exhibit precursor effects, all of which are 
found to be in good agreement ( Larkin and Varl amovl 1200 ih 
with fluctuation theory. As pseudogap effects become more 
pronounced with underdoping much of the fluctuation behav- 
ior appears to set in at a higher temper ature scale associated 
with T*, but oft en some fraction thereof ( iTan and Levinll2004l 
IXuef aal2000l) . 

Figures [O] and make the important point that precur- 
sor effects in the transverse thermoelectric response (axy) and 
conductivity a appear at higher temperatures (cx T*) as pseu- 
dogap effects become progressively more important; the dot- 
ted lines which have the strongest pseudogap continue to the 



highest temperatures on both Figures. Moreover, both trans- 
port coefficients evolve smoothly from a regime where they 
are pr esumed to be described by conventional fl uctuations 
(Larki n and Varlamovll200lllUssishkin et a/.Ll2002h as shown 
by the solid lines in the figures into a regime where their be- 
havior is ass ociated with a pseudogap. A n umber of people 
have argued dCorson ef 'an il999llXu et ali 2000) that nor- 
mal state vortices are responsible for the so called anomalous 
transport behavior of the pseudogap regime. We will argue in 
this Review that these figures may alternatively be interpreted 
as suggesting that bosonic degrees of freedom, not vortices, 
are responsible for these anomalous transport effects. 



C. Introduction to High Tc Superconductivity: Mott 
Physics and Possible Ordered States 

Most workers in the field of high Tc superconductivity 
would agree that we have made enormous progress in char- 
acterizing these materials and in identifying key theoretical 
questions and constructs. Experiment al progress, in large 
part, comes front transp ort studies (Tallon and LoraniLl200lL 
iTimusk and Stattlfi999h in addition to three powe rful spectro- 
scopies: photoemission fCampuzano ef fl/.,'2004; Damascelfi 
et ai, 2003), neutron ( AepDli et ai, 1997a b; Cheong et all 
[l991; Fong et a d Il99jfc IK astner et ai, 1998; Mook et a/J, 
1998; Rossat-Mignodef a/., 1991: Tranauada efaL. 1992.) and 
Josephson interferometry (Mathai et ai, 19951: iTsuei et al.{ 
1994; Woll man et al.l\199^ . These data have provided us 
with important clues to address related theoretical challenges. 
Among the outstanding theoretical issues in the cuprates are 
(i) understanding the attractive "mechanism" that binds elec- 
trons into Cooper pairs, (ii) understanding the evolution of the 
normal phase from Fermi liquid (in the "overdoped" regime) 
to marginal Fermi liquid (at "optimal" doping) to the pseudo- 
gap state, which is presumed to occur as doping concentration 
X decreases, and (iii) understanding the nature of that super- 
conducting phase which evolves from each of these three nor- 
mal states. 

The theoretical community has concentrated rather exten- 
sively on special regions and x-dependences in the phase dia- 
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FIG. 13 Transverse thermoelectric response (which relates to the 
Nernst coefficient) plotted here in fluctuation regime above Tc. Here 
OP and UD correspond to optimal and underdoping, respectively. 
Data were taken from.Xu et a/... ,2000. 
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gram which are presumed to be controlled by "Mott physics". 
There are different viewpoints of precisely what constitutes 
Mott physics in the metallic phases, and, for example, whether 
or not this is necessarily associated with spin-charge separa- 
tion. For the most part, one presumes here that the insulating 
phase of the parent compound introduces strong Coulomb cor- 
relations into the doped metallic states which may or may not 
also be associated with strong antiferromagnetic correlations. 
There is a recent, rather complete review (Lee et ai, 200|) 
of Mott physics within the scenario of spin-charge separation, 
which clarifies the physics it entails. 

Concrete experimental signatures of Mott physics are the 
observations that the superfluid density ps{T = 0, x) ^ 
as x 0, as if it were reflecting an order parameter for a 
metal insulator transition. More precisely it is deduced that 
Ps(0, x) oc X, at low X. Unusual effects associated with 
this linear- in- x dependence also show up in other experi- 
ments, such as the weight of coherence features in photoemis- 
sion data fCampuzano et ai, 2004; Damascelli et ai, 2003), 
as we ll as in thermodynamical signatures (Tallon and Loram. 
1200 Ih . 

While there is no single line of reasoning associated with 
these Mott constr aints, the low value of the su perfluid density 
has been argued dEmerv and Kivelsonlll995l) to be responsi- 
ble for soft phase fluctuations of the order parameter, which 
may be an important contributor to the pseudogap. However, 
recent concerns about this "phase fluctuatio n scenario" for the 
origin of the pseudogap have been r aised (Lei l2004 . It is 
now presumed by a number of groups that phase fluctuations 
alone may not be adequate and an additional static or fluctu- 
ating order of one form or another needs to be incorporated. 
Related to a competing o r co-existing order are conjectures 
(IChakravartv et ali 1200 i h that the disappearance of pseudo- 
gap effects around x ~ 0.2 is an indication of a "quantum 
critical point" associated with a hidden order parameter which 
may be responsible for the pseudogap. Others have associated 
small (Voita et ai, 2000) x or alternatively optimal (Varma, 
[^99) X with quantum critical points of a different origin. The 
nature of the other competing or fluctuating order has been 



conjectured to b e RVB llAnders on et all '(2004*), "d- density 
wave" ( Chakrav artv et g/.F 200 ll) . stripes ( Castellani et all 
ll995l:IZacha r et ai, 1997) or possib ly antiferromag netic spin 
flucttiations (tChubukov et aZ...1996i:lDemler and Zhang:.. 1 9951 
|PinesLll997h . The latter is another residue of the insulating 
phase. 

What is known about the "pairing mechanism"? Some 
would argue that this is an ill-defined question, and that su- 
perconductivity has to be understood through condensation 
energy arguments based for example on the data generated 
from the extrapolated normal state entropy dLoram et all 
11994 iTallon an d Loramll200ll) discussed above in the con- 
text of Fig. ^2 These condensation arguments are often built 
around the non-trivial assumption that there is a thermody- 
namically well behaved but meta-stable normal phase which 
coexists with the superconductivity. Others would argue that 
Coulomb effects are responsible for d-wave pairing, either di- 
rectly ( Leggett. 1999; Liu and Levin... 1997 , ). or indirectly vi a 
magnetic fluctuations (|Pineslll997l ; IScalapino et all Il98d) . 
Moreover, the extent to which the magnetism is presumed 
to persist into the metallic phase near optimal doping is con- 
troversial. Initially, NMR mea surements were interpreted as 
suggesting dMiUis ef aZ.Lll990l) strong antiferromagnetic fluc- 
tuations, while neutron measurements, which are generally 
viewed as the more conc lusive, do not provide compelling 
evidence dSi et all 1 19931) for their presence in the normal 
phase. Nevertheless, there are interesting neutron-me asured 
magnetic signatures (Kao et ai. 2000; Liu et ai. 1995; Zha 
et g/.. 1993) below Tc associated with rf-wave superconductiv- 
ity. There are also analogous STM effects which are currently 
of interest dWang and Leeil2003l) . 



D. Overview of the BCS-BEC Picture of High 
Superconductors in the Underdoped Regime 

We summarize here the physical picture associated with 
the more controversial BCS-BEC crossover approach to high 
Tc superconductors, based on the standard ground state of 
Eq. (0. One presumes that the normal state (pseudo)gap 
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arises from pre-formed pairs. This gap A(T) is thus present 
at the onset of superconductivity, and reflects the fact that 
Cooper pairs form before they Bose condense, as a result of 
sufficiently strong attraction. Below Tc, A evolves smoothly 
into the excitation gap of the superconducting phase. 

This behavior reflects the fact that Tc corresponds to the 
condensation of the q = pairs, while finite momentum pairs 
persist into the order phase as excitations of the condensate. 
As the temperature is progressively lowered below Tc, there 
are fewer and fewer such bosonic excitations. Their number 
vanishes precisely at T = 0, at which point the system reaches 
the ground state of Eq. Q. Here the superconducting or- 
der parameter and the excitation gap become precisely equal. 
This approach is based on a mean field scheme and is thus not 
equivalent to treating phase fluctuations of the order parame- 
ter. In general, below Tc, the superconducting order parameter 
Asc and the pseudogap contributions add in quadrature in the 
dispersion of Eq. (|6j as 

A\T) ^ Al{T) + Al^{T) (10) 

where Apg ^ as T ^ 0; this behavior should be contrasted 
with that in Eq. Q. 

Transport and thermodynamics in the normal and supercon- 
ducting phases are determined by two types of quasiparticles: 
fermionic excitations which exhibit a gap A(T), and noncon- 
densed pairs or bosonic excitations. Contributions from the 
former are rather similar to their (low T) counterparts in BCS 
theory, except that they are present well above Tc as well, 
with substantially larger A than in strict BCS theory. In this 
way, experiments which probe the opening of an excitation 
gap (such as Pauli susceptibility in the spin channel, and re- 
sistivity in the charge channel) will exhibit precursor effects, 
that is, the gradual onset of superconductivity from around T* 
to below Tc- 

The contributions to transport and thermodynamics of the 
bosonic excitations are rather similar to those deduced from 
the time dependent Ginsburg-Landau theory of superconduct- 
ing fluctuations, provided that one properly extends this the- 
ory, in a quantum fashion, away from the critical regime. It 
is these contributions which often exhibit divergences at Tc, 
but unlike conventional superconductors, bosonic contribu- 
tions (pre-formed pairs) are present, up to high temperatures 
of the order of T* /2, providing the boson density is suffi- 
ciently high. Again, these bosons do not represent fluctua- 
tions of the order parameter, but they obey a generically sim- 
ilar dynamical equation of motion, and, in the more general 
case which goes beyond Eq. Q, they will be coupled to order 
parameter fluctuations. Here, in an oversimplified sense there 
is a form of spin-charge separation. The bosons are spin-less 
(singlet) with charge 2e, while the fermions carry both spin 
and charge, but are non-Fermi-liquid like, as a consequence 
of their energy gap. 

In the years following the discovery of high Tc supercon- 
ductors, attention focused on the so-called marginal Fermi 
liquid behavior of the cuprates. This phase is associated 
with near-optimal doping. The most salient signature of this 
marginality is the dramatically linear in temperature depen- 
dence of the resistivity from just above Tc to very high tem- 



peratures, well above room temperature. Today, we know that 
pseudogap effects are present at optimal doping, so that the 
precise boundary between the marginal and pseudogap phases 
is not clear. As a consequence, addressing this normal state 
resistivity in the BCS-BEC crossover approach has not been 
done in detail, although it is unlikely to yield a strictly linear T 
dependence at a particular doping. Nevertheless, the bosonic 
contributions lead to a decrease in resistivity as temperature 
is lowered, while the opening of a d-wave gap associated 
with the fermionic terms leads to an increase, which tends to 
compensate. These are somewhat analogous to Aslamazov- 
Larkin and Maki-Thompson fluctuation contributions (Larkin 
and Varlamov. 2001) . The bosons generally dominate in the 
vicinity of Tc, but at higher T, in the vicinity of T*, the gapped 
fermions will be the more important. 

The vanishing of Tc with sufficiently large T* ( extreme 
underdoping) occurs in an interesting way in this crossover 
theory. This effect is associated with the localization of the 
bosonic degrees of freedom. Pauli principle effects in con- 
junction with the extended size of d-wave pairs inhibit pair 
hopping, and thus destroy superconductivity. One can specu- 
late that, just as in conventional Bose systems, once the bosons 
are localized, they will exhibit an alternative form of long 
range order. Indeed, "pair density wave" models have been 
invoked (Chen et ai. 2004) to exp lain recen t STM data (Fu 
et ai, 2004: IVershinin et al.l l2004t) . In contrast to the cold 
atom systems, for the d-wave case the BEC regime is never 
accessed. This superconductor-insulator transition takes place 
in the fermionic regime where /i « O.SEp. 

In the context of BCS-BEC crossover physics, it is not es- 
sential to establish the source of the attractive interaction. For 
the most part we skirt this issue in this Review by taking the 
pairing onset temperature T*{x) as a phenomenological in- 
put parameter. At the simplest level one may argue that as the 
system approaches the Mott insulating limit (x — > 0) fermions 
become less mobile; this serves to increase the effectiveness 
of the attractive interaction at low doping levels. 

The variable x in the high Tc problem should be viewed 
then, as a counterpart to the magnetic field variable in the 
cold atom problem which tunes the system through a Fesh- 
bach resonance; both x and field strength modulate the size 
of the attraction. However, for the cuprates, in contrast to the 
cold atom systems, one never reaches the BEC limit of the 
crossover theory. 

It is reasonable to presume based on the evidence to date, 
that pairing ultimately derives from Coulombic effects, not 
phonons, which are associated with / ~ pairing. While 
the widely used Hubbard model ignores these effects, longer 
ranged screened Coulomb interactions have been found (Liu 
and Levin, 12^997 ) to be attractive for electrons in a d-wave 
channel. In this context, it is useful to note that, similarly, 
in Hc^ short range repulsion destroys s-wave pairing, but 
leads to attraction in a hiaher (/ = 1) channel (Anderson 
and Brinkman, 1973). There is, however, no indication of 
pseudogap phenomena in He^, so that an Eliashberg extended 
form of BCS theory appears to be adequate. Eliashberg the- 
ory is a very different form of "strong coupling" theory from 
crossover physics, which treats in detail the dynamics of the 
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mediating boson. Interestingly, there is an upper bound to Tc 
in both schemes. For Eliashberg th eory this arises from th e 
induced effective mass corrections llLevin and VallsL Il983h . 
whereas in the crossover problem this occurs because of the 
presence of a pseudogap at Tc- 

Other analogies have also been invoked in comparing the 
cuprates with superfluid He'^, based on the propensity towards 
magnetism. In the case of the quantum liquid, however, there 
is a positive correlation betw een short range magnetic order- 
ing (i Levin and Vallsill983h and superconductivity while the 
correlation appears negative for the cuprates ( and some of 
their heavy fermion "cousins") where magnetism and super- 
conductivity appear to compete. 
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FIG. 15 A schematic representation of the Feshbach resonance for 
"Li. The effective interaction between the scattering of atoms in the 
triplet state is greatly enhanced, when the energy level of the singlet 
bound state approaches the continuum of the S = 1 channel. 



1 . Searching for the Definitive High Tc Theory 

In no sense should one argue that the BCS-BEC crossover 
approach precludes consideration of Mott physics. Nor should 
it be viewed as superior. Indeed, many of the experiments we 
will address in this Review have alternatively been addressed 
within the Mott scenario. Presumably both components are 
important in any ultimate theory. We have chosen in this Re- 
view to focus on one component only, and it is hoped that a 
comparably detailed review of Mott physics will be forthcom- 
ing. In this way one could compare th e strengths and weak- 
nesses of both schemes. A summary (Carlson et al. 
of some aspects of the Mott picture with emphasis on one di- 
mensional stripe states should be referred to as a very useful 
starting point. Also available now is a new Review fL ee et aUL 
12004) of Mott physics, based on spin charge separation. 

Besides those discussed in this Review, there is a litany of 
other experiments which are viewed as requiring explanation 
in order for a high Tc theory to be taken seriously. Among 
these are the important linear resistivity in the normal state 
(at optimal doping) and unusual temperature and x depen- 
dences in Hall data. While there has been little insight offered 
from the crossover scheme to explain the former, aspects of 
the latter have been succ essfully addresse d within pre-formed 
pair scenarios (Geshken bein et al 1 1199-1) . Using fermiology 
based schemes the major features of Dhotoemission(Norman 
and Pepin, 2003) and of neutron scattering and NM R data 
were studied theoretically and rather early on ( Ka o et al\ 
l2000t ISi et ai[ I1993L Eia et al\. Il993). along with a focus 
on the now famous 41 meV feature (Liu et al., 1995). These 
and other related approaches to the magnetic data, make use 
of a microscopically derived dSi et al\ Il992 l) random phase 
approximation scheme which incorporates the details of the 
Fermi surface shape and d-wave gap. Among remaining ex- 
periments in thi s litany is the origin o f incoherent transport 
along the c-axis dRoio and Levinlll993l) . 

On the basis of the above cited literature which addressed 
magnetic and c-axis data we are led to conclude that pseu- 



Here the afe-plane refers to the crystal plane in parallel with the layers of the 
cuprate compounds, and the c-axis refers to the axis normal to the layers. 



dogap effects and related crossover physics are not central to 
their understanding. The presence of a normal state gap will 
lead to some degree of precursor behavior above Tc, but d- 
wave nestin g relat ed effects such as reflected in neutron data 
dKao et ad 12000^. 'Li u et aIl ll995VZha ef a/.L Il993h wifl be 
greatly weakened a bove Tc by the blurring out of the d-wav e 
gap nodes ( Campuzano et a/.Ll2004HDamascelli et al I I2OOI . 
It should be said, however, that c-axis optical data do appear 
to reflect pseudogap effects (loffe and Millis, 1999) above Tc- 
(The a6-plane counterparts are discussed in Section lyLCI l. In 
summary, we are far from the definitive high Tc theory and the 
field has much to gain by establishing a collaboration with the 
cold atom community where BCS-BEC crossover as well as 
Mott physics ( Hofstetter et a/.. 2002) can be further elucidated 
and explored. 



E. Many Body Hamiltonian and Two Body Scattering 
Theory: Mostly Cold Atoms 

The existence of Feshbach resonances in a gas of ultracold 
fermionic atoms allows the interaction between atoms to be 
tuned arbitrarily, via application of a magnetic field. Figure 
^] presents a characteristic plot of the two body interaction 
potentials which are responsible for this resonance. These 
interaction potentials represent an effective description of all 
Coulomb interactions between the electrons and nuclei. For 
definiteness, we consider the case of ^Li here. In the region of 
magnetic fields near the "wide" Feshbach resonance, w 837G, 
we may focus on the two lowest atomic levels characterized 
by |F, tof), where the two quantum numbers correspond to 
the total spin and its z component. We specify the two states 
|1) = \\,\) and |2) = |i, —\)- In place of these quantum 
numbers one can alternatively characterize the two states in 
terms of the projections of the electronic (m^) and nuclear 
{rrii) spins along the magnetic field. In this way the states 1 
and 2 have nig = — i and = 1 and 0, respectively. 

The curve in Fig. [21 labelled 5 = 1 represents the inter- 
action between these two states. This coiTesponds to the so- 
called "open" channel and is associated with a triplet (elec- 
tronic) spin state. A Feshbach resonance requires that there 
be an otherwise unspecified near-by state which we refer to as 
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in the "closed" channel. The potential which gives rise to this 
singlet (electronic) spin state is shown by the S — curve in 
Fig. ^] A Feshbach resonance arises when this bound state 
lies near the zero of energy corresponding to the continuum of 
scattering states in the triplet or open channel. The electronic 
Zeeman coupling to the triplet state makes is possible to shift 
the relative position of the two potentials by application of a 
magnetic field. This shift is indicated by 5 in Fig. [21 In this 
way, the resonance is magnetically tunable and the scattering 
length of the system will vary continuously with field, diverg- 
ing at that magnetic field which corresponds to the Feshbach 
resonance. We represent the deviation of this field from the 
resonance condition by the "detuning parameter" i^. When i/ 
is large and positive, the singlet state is substantially raised 
relative to the triplet- and it plays relatively little role. By 
contrast when the magnitude of ly is large, but with negative 
sign, then the physics is dominated by the singlet or closed 
channel. 

The most general form for the Hamiltonian contains two 
types of particles: fermions which represent the open chan- 
nel and "molecular bosons" which represent the two fermion 
bound state of the closed channel. We will also refer to the 
latter as "Feshbach bosons". These, in turn, lead to two types 
of interaction effects: those associated with the direct attrac- 
tion between fermions, parameterized by U, and those asso- 
ciated with "fermion-boson" interactions, whose strength is 
governed by g. The latter may be viewed as deriving from hy- 
perfine interactions which couple the (total) nuclear and elec- 
tronic spins. In the high Tc applications, one generally, but 
not always, ignores the bosonic degrees of freedom. Other- 
wise the Hamiltonians are the same, and can be written as 

k,cr q 

+ 5Z ^('''^')"q/2+k,T'^q/2-ka'^q/2-k'a"q/2+k\T 
q,k,k' 

+ (■9(k)fe]j«q/2-k,i"q/2+k.T + ^■'') (11) 

q,k 

Here the fermion and boson kinetic energies are given by 
ej(. — fc2/(2m), and e™ = q^/(2M), and ly is an impor- 
tant parameter which represents the magnetic field "detun- 
ing". M — 2m is the bare mass of the Feshbach bosons. 
For the cold atom problem a —1,1 represents the two dif- 
ferent hyperfine states_£Ll}_^n^j2}}^_Onejna^ignore (Tim- 
mermans et ai, 2001 ) the non-degeneracy of these two states, 
since there is no known way for state 2 to decay into state 1 . 
In the two channel problem the ground state wavefunction is 
slightly modified and given by 

'i'o = *o ® (12) 
where the normalized molecular boson contribution '^'q is 

as discussed originally by iKostvrko and RanningeA Il996l 
Here A is a variational parameter We present the resulting 
variational conditions for the ground state in AppendixIXl 
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FIG. 16 Characteristic behavior of the scattering length as in the 
three regimes, as diverges when i^o — 2fi. Here, and in subsequent 
figures ff) is always in units of Ef and we estimate 1 G ~ 60 Ef for 
a typical Ef = 2 ^K, for either ''Li or "'"K. 



Whether both forms of interactions are needed in either 
the high Tc or cold atom systems is still under debate. The 
bosons (5y) of the cold atom problem ('Holl and et all 1200 it 
Timmermans ef gii, 2001) will be referred to as Feshbach 
bosons. These represent a separate species, not to be con- 
fused with the fermion pair (al_a^_iJ operators. Thus we call 
this a "two channel" model. In this review we will discuss 
the behavior of crossover physics both with and without these 
Feshbach bosons (FB). Previous studies of high Tc supercon- 
ductors have i nvoked a sim ilar bosonic ter m (Geshkenbein 
et a/.. ll997l:lMicnas et adll990l:lRanninger and Robirifl99l 
as well, although less is known about its microscopic ori- 
gin. This fermion-boson coupling is not to be confused with 
the coupling between fermions and a "pairing-mechanism"- 
related boson ([6-1- b^a^a) such as phonons. The coupling 
b^aa and its conjugate represents a form of sink and source for 
creating fermion pairs, inducing superconductivity in some 
ways, as a by-product of Bose condensation. We will em- 
phasize throughout that crossover theory proceeds in a similar 
fashion with and without Feshbach bosons. 

It is useful at this stage to introduce the s-wave scattering 
length, tts- This parameter is associated with the scattering 
amplitude in the low energy limit, for a two body scattering 
process. This scattering length can, in principle, be deduced 
from Figure \^ It follows that is negative when there is 
no bound state, it tends to — oo at the onset of the bound state 
and to +00 just as the bound state stabilizes. It remains posi- 
tive but decreases in value as the interaction becomes increas- 
ingly strong. The magnitude of is arbitrarily small in both 
the extreme BEC and BCS limits, but with opposite sign (See 
Fig. ll6l l. The fundamental postulate of crossover theory is that 
even though the two-body scattering length changes abruptly 
at the unitary scattering condition (\as\ = oo), in the N-body 
problem the superconductivity varies smoothly through this 
point. Throughout this Review, we presume an equivalence 
between the "pseudogap" and unitary scattering regimes. 

An important goal is to understand how to incorporate all 
the constraints of the two body physics into the superfluid gap 
equation. As a by-product, one avoids implicit divergences 
which would also appear in the problem. To begin we ignore 
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Feshbach effects. The way in which the two body interaction 
U enters to characterize the scattering (in vacuum) is different 
from the way in which it enters to characterize the N-body 
processes leading to superfluidity. In each case, however, one 
uses a T-matrix formulation to sum an appropriately selected 
but infinite series of terms in U. For the two-body problem in 
vacuum, we introduce the scattering length, Og, 



TO 



1 



(14) 



Aira 

which is related to U via the Lippmann-Schwinger equation 

m _ 1 

" U 



47ra, 



E- 



2ek 



(15) 



In this way one may solve for the unknown U in terms of 
fls or Uq. The two-body T-matrix equation [Eq. (I15H can be 
rewritten as 



u = TU„, r = ( 1 + ^ 



where we define the quantity Uc as 



2ek 



(16) 



(17) 



Here Uc is the critical value of the potential associated with the 
binding of a two particle state in vacuum. Specific evaluation 
of Uc requires that there be a cut-off imposed on the above 
summation, associated with the range of the potential. 

Now let us turn to the Feshbach problem, where the super- 
fluidity is constrained by a more extended set of parameters 
U, g, V and also /i. Provided we redefine the appropriate "two 
body" scattering length. E quation J15t holds even in the pres- 
ence of Feshbach effects ("Milstein ef al). '2002': Ohashi and 
Griffin, 2002). It has been shown that U in the above equa- 
tions is replaced by C/g/ / and by a* 
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(18) 



Here it is important to stress that many body effects, via the 
fermionic chemical potential ^ enter into a*. 

More precisely the effective interaction between two 
fermions is Q dependent. It arises from a second order pro- 
cess involving emission and absorption of a molecular boson. 
The net effect of the direct plus indirect interactions is given 
by 



UcffiQ) = U + g^DoiQ) , 



where 



DoiQ) = i/[in. 



(19) 



(20) 



is the non-interacting molecular boson propagator What ap- 
pears in the gap equation, however, is 



UcffiQ^O)^Ucff^U 



(21) 



Experimentally, the two body scattering length a* varies 
with magnetic fi eld B. Near re sonance, it can be parameter- 
ized ( Kokkelma ns et aZ.L l2002h in terms of three parameters 
of the two body problem Uq, go and the magnetic field which 
at resonance is given by Bq. One fits the scattering length to 
af,g(l — g '^g ) so that the width w is related to go and the 
background scattering length ai,g is related to Uq. Away from 
resonance, the crossover picture requires that the scattering 
length vanish asymptotically as 



0, 



in extreme BEC limit. 



(22) 



so that the system represents an ideal Bose gas. Here it should 
be noted that this BEC limit corresponds to ly — > The 
counterpart of Eq. (I22> also holds for the one channel prob- 
lem. Note that the actual resonance of the many body problem 
will be somewhat shifted from the atomic-fitted resonance pa- 
rameter Bq. 

It is convenient to define the parameter i/q which is directly 
related to the difference in the applied magnetic field B and 
Bo 



(23) 



where A/i° is the difference in the magnetic moment of the 
singlet and triplet paired states. We will use the parameter vq 
as a measure of magnetic field strength throughout this Re- 
view. The counterpart of Eq. MA\ is then 
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(24) 



where Uq 



and 5q = Uo^fJ.'^w. For the purposes of 



this Review we will use the scattering length a*, although for 
simplicity we often drop the asterisk and write it as a^. 

We can again inv ert the Lippmann-Schwinge r or T-matrix 
scattering equation ("Bruun and Pethicl^, '(2004'; Kokkelmans 
et a/.'T2002) to arrive at the appropriate parameters g and ly 
which enter into the Hamiltonian and thus into the gap equa- 
tion. Just as in Eq. (O, one finds a connection between mea- 
surable parameters (with subscript 0) and their counterparts in 
the Hamiltonian (without the subscript): 



U = TUo, g^Tgo, v - vq ^ -T 



Uc 



(25) 



To connect the various energy scales, which appear in the 
problem, typically IG « GOEp for both ^Li and ^Or. Here 
it is assumed that the Fermi energy is Ep w 2 /xK. 

In Fig. ll6l and using Eq. M%\ . we plot a prototypical scat- 
tering length kptts = kpa* as a function of the magnetic field 
dependent parameter, vq. This figure indicates the BEC, BCS 
and PG regimes. Here, as earlier, the PG regime is bounded 
on the left hy fi — Q and on the right by ^{Tc) ~ 0. 

A key finding associated with this plot is that the PG regime 
begins on the so-called "BEC side of resonance ", indepen- 
dent of whether Feshbach effects are included or not. That is, 
the fermionic chemical potential reaches zero while the scat- 
tering length is positive. This generic effect of the ground 
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state self consistent equations arises from the Pauli principle. 
For an isolated system, two particle resonant scattering occurs 
at U = Uc, where Uc is defined in Eq. (I17t . In the many 
body context, because of Pauli principle repulsion between 
fermions, it requires an effectively stronger attraction to bind 
particles into molecules. Stated alternatively, the onset of the 
fermionic regime (/i > 0) occurs for more strongly attractive 
interactions, than those required for two body resonant scat- 
tering. As a consequence of /i > 0, there is very little conden- 
sate weight in the Feshbach boson channel near the unitary 
limit, for values of the coupling parameter g^, appropriate to 
currently studied (and rather wide) Feshbach resonances. 

Finally, it is useful to compare the two interaction terms 
which we have defined above, (J7e// and [/*), and to con- 
trast their behavior in the BEC and unitary limits. The quan- 
tity U* is proportional to the scattering length; it reflects the 
two body physics and necessarily diverges at unitarity, where 
Ueff = Uc- By contrast, in the BEC regime, U*, or a* ap- 
proaches zero, (under the usual presumption that the interac- 
tion is a contact interaction). This vanishing of U* should be 
compared with the observation that the quantity Ueff, which 
reflects the many body physics, must diverge in this BEC 
limit. 



1 . Differences Between One and Two Channel Models: 
Physics of Feshbach bosons 

The main effects associated with including the explicit Fes- 
hbach resonance are twofold. These Feshbach bosons provide 
a physical mechanism for tuning the scattering length to be 
arbitrarily large, and in the BEC limit, they lead to different 
physics from the one channel problem. As will be illustrated 
later via a comparison of Figures |23] and |25] in an artificial 
way one can capture most of the salient physics of the uni- 
tary regime via a one channel approach in which one drops all 
molecular boson related terms in the Hamiltonian, but takes 
the interaction U to be arbitrary. In this way there is very little 
difference between fermion-only and fermion-boson models. 
There is an important proviso, though, that one is not dealing 
with narrow resonances which behave rather differently near 
unitarity, as will be seen in Fig.|32] In addition, as in the one 
channel model a form of universality (Ho, 2004) is found in 
the ground state at the unitary limit, provided g is sufficiently 
large. This is discussed later in the context of Fig. 1301 

We emphasize that, within the BEC regime there are three 
important effects associated with the Feshbach coupling g. As 
will become clear later, (i) in the extreme BEC limit when g 
is nonzero, there are no occupied fermionic states. The num- 
ber constraint can be satisfied entirely in terms of Feshbach 
bosons, (ii) The absence of fermions will, moreover, greatly 
weaken the inter-boson interactions which are presumed to 
be mediated by the fermions. In addition, as one decreases 
|i7e//| from very attractive to moderately attractive (i.e., in- 
creases a* on the BEC side) the nature of the condensed pairs 
changes. Even in the absence of Feshbach bosons (FB), the 
size of the pairs increases. But in their presence the admix- 
ture of bosonic (i.e., molecular bosons) and fermionic (i.e.. 



Cooper pair) components in the condensate is continuously 
varied from fully bosonic to fully fermionic. Finally (iii) the 
role of the condensate enters in two very different ways into 
the self consistent gap and number equations, depending on 
whether or not there are FB. The molecular boson condensate 
appears principally in the number equation, while the Cooper 
pair condensate appears principally in the gap equation. To 
make this point clearer we can refer ahead to Eq. ( I81> which 
contains the FB condensate contribution By contrast the 
Cooper condensate contribution Age enters into the total ex- 
citation gap A as in Eq. dlOl l and A is, in turn, constrained 
via Eq. \69\ . Points (i) and (iii), which may appear some- 
what technical, have essential physical implications, among 
the most important of these is point (ii) above. 

We stress, however, that /or the PG and BCS regimes the 
differences with and without FB are, however, considerably 
less pronounced. 



F. Current Summary of Cold Atom Experiments: 
Crossover in the Presence of Feshbach Resonances 

There has been an exciting string of developments over the 
past few years in studies of ultracold fermionic atoms, in par- 
ticular, ^Li and ^"K, which have been trapped and cooled via 
magnetic and optical means. Typically these traps contain 
^ 10^ atoms at very low densities w lO^'^ cm"'^ when cooled 
to very low T. Here the Fermi temperature Ep in a trap can be 
estimated to be of the order of a mi cro-Kelvin. It was a rgued 
on the basis of BCS theory alone ( Houb iers et a/.L [1997 1. and 
rather early on (1997), that the temperatures associated with 
the superfluid phases may be attainable in these trapped gases. 
This set off a search for the "holy grail" of fermionic super- 
fluidity. That a Fermi degenerate state could be reached at 
all is itself quite remarkable; this was was first reported (De- 
Mai-co and Jin. 1999 ) by Jin and deMarco in 1999. By late 
2002 reports of unusual hydrodynamics in a degenerate Fermi 
gas indicated that strong interactions were present (O'Hara 
et ai, 2002). This strongly interacting Fermi gas (associated 
with the unitary scattering regime) has attracted widespread 
attention independent of the search for superfluidity, because 
it appears to be a prototype for analogous systems in nuclear 
physics (Baker, 1999; Heiselberg, 2004a) and in quark-gluon 
plasmas ( It akura^ .2003; .Kolb and Heinz, .2003>) . Moreover, 
there has been a fairly extensive body o f analytic work on the 
ground state properties of this regime (' Carlson et allll999t 
■Heiselberg. 2001..) . which goes beyond the simple mean field 
wave function ansatz. 

As a consequence of attractive s-wave interactions between 
fermionic atoms in different hyperfine states, it was antici- 
pated that dimers could al so be made. Indeed, these molecules 
formed rather efficiently ("Cubizolles et a/.'. '2003': 'Regal et all 
2003; Strecker et ai, 2003) as reported in mid-2003 either via 
three body recombination flJochim et a/.Ll2003al) or by sweep- 
ing the magnetic field across a Feshbach re sonance. More- 
over, they are extremely long lived (Streck er et al.l E003h . 
From this work it was relatively straightforward to anticipate 
that a Bose condensate would also be achieved. Credit goes 
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FIG. 17 (Color) Spatial density profiles of a molecular cloud of 
trapped ''"K atoms in the BEC regi me in the transverse directions 
after 20 ms of free expansion rfrom lOreiner et a i. 2003), showing 
thermal molecular cloud above Tc (left) and a molecular condensate 
(right) below Tc. (a) shows the surface plots, and (b) shows the cross- 
sections through images (dots) with bimodal fits (lines). 

to theorist s such as Holland (|Holland et a/.Ll200ll) and to Tim- 
mermans dTimmermans ef am20o1F and their co-workers for 
recognizing that the superfluidity need not be only associated 
with condensation of long lived bosons, but in fact could also 
derive, as in BCS, from fermion pairs. In this way, it was 
argued that a suitable tuning of the attractive interaction via 
Feshbach resonance effects would lead to a realization of a 
BCS-BEC crossover 

By late 2003 to early 2004, four group s teourd el et al\. 
l2004 iGreiner et a/.L l2003HTochim et aZ.L r2 003b: Zwierlein 
et ai, 2003) had observed the "condensation of weakly bound 
molecules" (that is, on the as > side of resonance), and 
shortly thereafter a number h ad also reported evidence for su- 
perfluidity on the BCS side (Chin et al. . '2004': ' Kinast et a/1 
|2004; Re gal et ai. 2004; Zwierl ein et ai. 2004j). The BEC 
side is the more straightforward since the presence of the su- 
perfluid is reflected in a bi-modal distribution in the density 
profile. This is shown in Fig. ll7l from (Greiner et ai, 2003), 
and is conceptually s imilar to the behavior for condensed Bose 
atoms l lDalfovo et a/.i.il999i) . On the BEC side but near reso- 
nance, the estimated Tc is of the order of 500 nK, with conden- 
sate fractions varying from 20% or so, to nearly 100%. The 
condensate lifetimes are relatively long in the vicinity of res- 
onance, and fall off rapidly as one goes deeper into the BEC. 
However, for < there is no clear expectation that the den- 
sity profile will provide a signature of the superfluid phase. 

These claims that superfluidity may have been achieved on 
the BCS side (og < 0) of resonance were viewed as partic- 
ularly exciting. The atomic community, for the most part, 
felt the previous counterpart observations on the BEC side 
were expected and not significantly different from conden- 
sation in Bose atoms. The evidence for this new form of 
"fermionic supe rfluidity" rests on studies dRegal et a/.Ll2004l 
IZwierlein et alx .2004.) that perform fast sweeps from nega- 



tive Us to positive across the resonance. The field sweeps 
allow, in principle, a pairwise projection of fermionic atoms 
(on the BCS side) onto molecules (on the BEC side). It is 
presumed that in this way one measures the momentum dis- 
tribution of fermionic atom pairs. The existence of a con- 
densate was thus inferred. Other experiments which sweep 
across the Feshbach resonance adiabatically, measure the size 
of the clou d after relea se (Bourdel et ai, 2004) or within a 
trap (Barte nstein et al\. ^004b). Evidence for superfluidity 
on the BCS side, which does not rely on the sweep experi- 
ments, has also been deduced from collective excitations of a 
fermionic gas (B artenstein et al. l l2004atlKiiiast et a/.Ll2004l) . 
Pairing gap measurements with radio frequency (RF) probes 
dChin et a/.L i2004) hav^imn^^beer^terDreted (Kinnunen 
et ai, 2004a) as evidence for superfluidity, although it appears 
more likely that these experiments establish only the existence 
of fermionic pairs. Quite recently, evidence for a phase tran- 
sition has been presented via thermodynamic measurements 
and accompanying theory (Kina st et ai, 2005). The latter, 
like the theorv (iKinnunen et a/.L l2004a) of RF experiments 
jChin et all l2004^ ■ is based on the formalism presented in 
this Review. 

Precisely what goes on during the sweep is not entirely un- 
derstood. It should be stressed again that where the scatter- 
ing length changes sign is not particularly important to the 
N-body physics. Thus starting on the "BCS side of reso- 
nance " and ending on the "BEC side of resonance " involves a 
very continuous sweep in which the physics is not qualitatively 
changed. It has been speculated that for this sweep procedure 
to work a large percentage of the Cooper pair partners must be 
closer than the inter-atomic spacing. Others have alternatively 
viewed the sweep as having little to do with the large Cooper 
pairs of BCS superconductivity. 

It is useful to establish the various pair sizes here. In the 
PG regime the normal state consists of a significant number 
of small (compared to BCS) pre-formed pairs. Below Tc, the 
condensate pair size is also small. Finally, there are excited 
pair states in the superfluid phase, with characteristic size ^pg. 
The Cooper pair size £ as p lotted in Fig.[^was deduced from 
the condensate (Randeria, 1995) at T = 0, but this can be 
shown to be rather similar to corresponding t o the size of 
excited pair states (IChenll200alChen e/anil998h bofli above 
and below Tc- It seems plausible that during the sweep there 
is some rearrangement of excited states (both fermionic and 
bosonic) and condensate. In addition, all pairs (excited and 
condensate) contract in size. When, at the end of a sweep, 
they are sufficiently small, then they are more "visible". What 
appears to be essential, however, is that the sweep time scales 
are rapid enough so that the existence of a condensate at the 
initial starting point is a necessary and sufficient condition for 
observing a condensate at the end of a sweep. 



G. T-Matrix-Based Approaches to BCS-BEC Crossover in 
the Absence of Feshbach Effects 

In this subsection we introduce the concept of a T-matrix 
within the pairing channel. This is defined in Eq. M6\ be- 
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and theories differ only on what is the nature of the pair 
susceptibility x(*9)^ ™d the associated self energy of the 
fermions. Here and throughout we use Q to denote a four- 
vector and write J^q = '^J2in X]q' where is a Mat- 
subara frequency. Below Tc one can also consider a T-matrix 
approach to describe the particles and pairs in the condensate. 
For the most part we will defer extensions to the broken sym- 
metry phase to Section HTHl 



1 . Review of BCS Theory Using tlie T-matrix Approacli 



FIG. 18 Characteristic behavior of Cooper pair size ^ in the con- 
densate in the three regimes. The dotted vertical lines are not sharp 
transitions but indicate where /i = (left line) and where A{Tc) ~ 
(right line). The dimensionless kpS, gives the ratio between pair size 
^ and the mean inter-particle spacing l/kp- Pairs shrink and become 
tightly bound in the deep BEC regime. 



low and it represents an effectively screened interaction in 
which the appropriate polarizability (x) is associated with the 
particle-particle rather than the particle-hole channel, as in the 
more conventional dielectric susceptibility. The T-matrix thus 
appears as an infinite series of two particle "ladder" diagrams. 

Away from zero temperature, a variety of different many 
body approaches have been invoked to address the physics of 
BCS-BEC crossover For the most part, these revolve around 
T-matrix schemes. Here one solves self consistently for the 
single fermion propagator (G) and the pair propagator (t). 
That one stops at this level without introducing higher order 
Green's functions (involving three, and four particles, etc) is 
believed to be adequate for addressing a leading order mean 
field theory such as that represented by Eq. Q. One can see 
that pair-pair (boson-boson) interactions are only treated in a 
mean field averaging procedure; they arise exclusively from 
the fermions and are presumed to be sufficiently weak so as 
not to lead to any incomplete condensation in the ground state, 
as is compatible with Eq. Q. 

One can view this approach as the first step beyond BCS in 
a hierarchy of mean field theories. In BCS, above Tc one in- 
cludes only the bare fermion propagator Gq. Below Tc, pairs 
play a role but only through the condensate. At the next level 
one accounts for the interaction between particles and noncon- 
densed pairs in both the normal and superconducting states. 
The pairs introduce a self energy E into the particles, which 
represents a correction to BCS theory. The pairs are treated 
at an effective mean field level. By truncating the equations 
of motion in this way, the effects of all higher order Green's 
functions are subsumed into t in an averaged way. 

Below we demonstrate that at this T-matrix level there are 
a number of distinct schemes which can be implemented to 
address BCS-BEC crossover physics. These same alternatives 
have also entered into a discussion of pre-formed pairs as they 
appear in treatments of superconducting fluctuations. Above 
Tc, quite generally one writes for the T-matrix 



u 



1 + UxiQ) 



(26) 



In order to assess alternative schemes, it is useful to re- 
view BCS theory within a T-matrix formalism. In BCS the- 
ory, pairs explicitly enter into the problem below Tc, but only 
through the condensate. These condensed pairs are associated 
with a T-matrix given by 

i,c(Q) - ^AlS{Q)/T (27) 

with fermionic self energy 

S,c (K) = (Q)^o - ^) (28) 

Q 

so that T.sc{K) = -A^^Go(-if). Here, and fliroughout, Gq 
is the Green's function of the non-interacting system. The 
number of fermions in a BCS superconductor is given by 



K 



and 



G{K) = [G^\K)-i:sc{K)]-' 



(29) 



(30) 



Doing the Matsubara summation in Eq. ( I29> . one can then 
deduce the usual BCS expression for the number of particles, 
which determines the fermionic chemical potential 



1 - 



fiEi.) 



where 



Eu 



(31) 



(32) 



We need, however, an additional equation to determine 
Asc{T). The BCS gap equation can be written as 



1 + t/XBCs(0) = 0, T<Tc 



(33) 



where 



XBCsiQ) - J2 G{K)G^{Q - K) (34) 



This suggests that one consider the uncondensed or normal 
state pair propagator to be of the form 



t{Q) 



u 



1 + UxBcsiQ) 



(35) 
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then in the superconducting state we have a BEC like condi- 
tion on the pair chemical potential /ipai,- defined by 

t^^{Q = 0) ~ fipair X const. (36) 

where the overall constant is unimportant for the present pur- 
poses. Thus we say that the pair chemical potential satisfies 

fJ-pair =0, T <Tc . (37) 

That fipair vanishes at all T < Tcis a stronger condition than 
the usual Thouless condition for Tc- Moreover, it should be 
stressed that BCS theory is associated with a particular asym- 
metric form for the pair susceptibility. These uncondensed 
pairs play virtually no role in BCS superconductors but the 
structure of this theory points to a particular choice for a pair 
susceptibility. We can then write the self consistent condition 
on Age which follows from Eq. ( I33t . Using Eqs. Ml\ I34> . 
we find the expected form for the BCS gap equation. This 
represents a constraint on the order parameter Ag^: 

A.,(T) = -f/^ A.,(T)i^||pil . (38) 

k 

The above discussion was presented in a somewhat differ- 
ent w av in a paper bv Kadanoff and Martin ( Kadanoff and 
Martin Il96ll) . 



2. Three Choices for the T-matrix of the Normal State: 
Problems with the Nozieres Schmitt-Rink Approach 

On general grounds we can say that there are three obvi- 
ous choices for x{Q) which appears in the general definition 
of the T-matrix in Eq. \26\ . Each of these has been studied 
rather extensively in the literature. All of these introduce cor- 
rections to BCS theory and all were motivated by attempts to 
extend the crossover ground state to finite T, or to understand 
widespread pseudogap effects in the high Tc superconductors. 
In analogy with Gaussian fluctuations, one can consider 

K 

with self energy 

^o{K) ^Y.^[Q)Go{Q - K) (40) 

Q 

which appears in G in the analogue of Eq. (I30> . The number 
equation is then deduced by using Eq. ( I29> . 

This scheme was adopted by Nozieres and Schmitt-Rink 
(NSR ) , alt hough these authors (Nozieres and Schmitt- Rin5 
Il985t jRanderia, 1995) approximated the number equation 
(ISerene', '1989) by using a leading order series for G in 
Eq. (EH) with 

G = Go + GoSoGo . (41) 



It is straightforward, however, to avoid thi s approximation in 
Dyson's eq uation, and a number of groups dJanko et a/.lll997|: 
IPieri et a7T E004) have extended NSR in this way. 
Similarly one can consider 

X{Q) ^Y.^^^)G^Q ~ ^) (42) 

K 

with self energy 

t{K) =Y,t{Q)G{Q ~ K) . (43) 

Q 

This latter scheme (sometimes known as fluctuation-exchange 
or FLEX) has been also extensively discussed in the lit- 
erature, by among ot hers, Ha ussmann (Haussmann, 1993j), 
Tchernyshy ov ( Tcher nvshvovt 1 1 9971) . Micnas a nd colleagues 
(iMicn as e7a/.i.il995i) and Yamada and Yanatse (lYanase et al\. 
2003). 

Finallv. we can contemplate the asvmmetric form (Chen 
et ai, 1 998.: Giovannini and Berthod.2001.) for the T-matrix, 
so that the coupled equations for t{Q) and G{K) are based on 

x(Q)=5]G(if)Go(Q-if) (44) 

K 

with self energy 

i:{K) =Y,t{Q)Go[Q - K) (45) 
Q 

Each of these three schemes determines the superfluid tran- 
sition temperature via the Thouless condition. Thus ^pair — 
leads to a slightly different expression for T^, based on the dif- 
ferences in the choice of T-matrix which appears in Eq. ( I36l l. 

A central problem with the NSR scheme (x ~ GqGq) is 
that it incorporates self energy effects only through the num- 
ber equation. The absence of self energy effects in the gap 
equation is equivalent to the statement that pseudogap effects 
only indirectly affect T^. the particles acquire a self energy 
from the pairs but these self energy effects are not fed back 
into the propagator for the pairs. On physical grounds one an- 
ticipates that the pseudogap must have a direct effect on Tc 
as a consequence of the associated reduction in the density of 
states. Oth er problems related to the thermodynamics were 
pointed out dSofo and Balseirol 'l992) when NSR was applied 
to a two dimensional system. Finally, because only the num- 
ber equation contains self energy effects, it is not clear if one 
can arrive at a proper vanishing of the superfluid density ps at 
Tc within this approach. This requires a precise cancellation 
between the diamagnetic current contributions (which depend 
on the number equation) and the paramagnetic terms (which 
reflect the gap equation). 

One might be inclined, then, to prefer (Seren4 Il989h the 
FLEX (x ~ GG) scheme since it is (/)-derivable, in the sense 
of Baym (Bavm, 1962). This means that it is possible to write 
down a closed form expression for the thermodynamical po- 
tential, from which one derives other physical quantities via 
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functional derivatives of this potential. Theoretical consis- 
tency issues in this appro ach have been rath er exhaustively 
discussed by Haussmann llHaussmanni Il993h above T^. We 
are not aware of a fully self consistent calculation of ps below 
Tc, at the same level of completeness as Haussmann's normal 
state analysis (or, for that matter, of the counterpart discussion 
to Section lTll. Al and accompanying AppendixlEt. However, in 
principle, it should be possible to arrive at a proper expression 
for the superfluid density, providi ng Ward identities are im- 
posedj^ier^s some ambiguitv (Benard et al.\'l993\ Micnas 
et a/. ni995n ^hernvshvov . 1997) about whether pseudogap 
effects are present in the FLEX approach; a consensus has not 
been reached at this time. Numerical work based on FLE X is 
more extensive (Benai'd et a/.L 1 1 9931 iDeisz et a/.L 1 1998h than 
for the other two alternative schemes. 

It will be made clear in what follows that, if one's goal is to 
extend the usual crossover ground state ofEq. to finite tem- 
peratures, then one must choose the asymmetric form for the 
pair susceptibility . Here x ~ GGq. This is different from the 
approach of Nozieres and Schmitt-Rink in which there is no 
evident consistency between their finite T treatment and the 
presumed ground state. Some have claimed ( Randeria. 19951) 
that an NSR-based approach, when extended below Tc, repro- 
duces Eqs. i69\ and ilOi . But elsewhere in the literature it 
has been shown that these ground state self consistency con- 
ditions are derive^Jroin^Jowe^JeveHheor^Kostyrko and 
Ranninger, 199^ They arise from a T = limit of what is 
called the saddle point approximation, or mean field approach. 
When pairing fluctuations bevond the saddle point (Nozieres 
and Schmitt-Rink, J 985) are included, the number equation is 
changed. Then th e ground state is no longer that of Eq. Q 
llPieri et a/ .'. 2004). In the cold atom literature, the distinction 
between the NSR-extrapolated ground state and the ground 
state of Eq. Q has been pointed ou t in the context of atomic 
density profiles llPerah ef alll2004bl) within traps. 

It should be evident from Eq. ( I34> and the surrounding 
discussion that the asymmetric form is uniquely associated 
with the BCS-like behavior implicit in Eq. Q- The other 
two T-matrix approaches lead to different ground states which 
should, however, be very interesting in their own right. These 
will need to be characterized in future. 



Additional support for the asymmetric form, which will be 
advocated here, was provided by Kadanoff and Martin. In 
their famous paper they noted that several people had sur- 
mised that the form for inv olving GG would be more ac - 
curate. However, as claimed in ( iKadanoff and Martirifl96ll) . 
"This surmise is not correct." Aside from theoretical counter- 
arguments which they present, the more symmetric combina- 
tion of Green's functions "can also be rejected experimentally 
since they give rise to a specific heat." Similar arguments in 
support of the asvmmetric form were introduced (Hassing and 
Wilkins. ll973l) into the literature in the context of addressing 
specific heat jumps. 



H. Superconducting Fluctuations: a type of Pre-formed 
Pairs 

While there are no indications of bosonic degrees of free- 
dom, (other than in the condensate), within strict BCS the- 
ory, it has been possible to access these bosons via probes 
of superconducting fluctuations. Quasi-one dimensional, or 
quasi-2D superconductors in the presence o f significant dis- 
order exhibit fluctuation effects ( fLarkin and Varlamovtl20()ll) 
or precursor pairing as seen in "paraconductivity", fluctuation 
diamagnetism, as well as other unusual behavior, often con- 
sisting of divergent contributions to transport. One frequently 
computes (Kwon and Dorsev. 1999) these bosonic contribu- 
tions to transport by use of a time dependent Ginzburg-Landau 
(TDGL) equation of motion. This is rather similar to a Gross- 
Pitaevskii (GP) formalism except that the "bosons" here are 
highly damped by the fermions. 

Alternatively T-matrix based approaches (involving all 
three choices of x{Q)) have been extensively used to discuss 
conventional superconducting fluctuations. The advantage of 
these latter schemes is that one can address both the anoma- 
lous bosonic and fermionic contributions to transport through 
the famous Aslamazov-Larki n and Maki-Thompson diagrams 
( iLarkin and Varlamovll200il) . We defer a discussion of these 
issues until AppendixiBl 

It is useful to demonstrate first how conventional supercon- 
ducting fluctuations behave at the lowest level of self con- 
sistency, called th e Hartree approxim ation. This scheme is 
closely associatedl'Staiic e t adl2003 ?) with a GGq T-matrix, 
as is discussed in more detail in Appendixicl It is also closely 
associated with BCS theory, for one can show that, in the spirit 
of Eq. Q, at this Hartree level the excitation gap A(Tc) and 
Tc lie on a specific BCS curve (specified by n and U). What is 
different from strict BCS theory is that the onset of supercon- 
ductivity takes place in the presence of a finite excitation gap 
(ie, pseudogap), just as shown in Fig. |3l This, then, reflects 
the fact that there are pre-formed pairs at Tc. By contrast with 
high Tc superconductors, however, in conventional fluctuation 
effects, the temperature T* at which pairs start to form is al- 
ways extremely close to their condensation temperature Tc. 
We thus say that there is a very narrow critical region. 

In Hartree approxirnatedGinzburg^andautheo^ 
and Wilkins. il973i) the free energy functional is given by 

F[^] = ao(T-T*)|*p+6|*|4 « ao(T-T*)|*p+26A2|*|2 

(46) 

Here plays the role of a pseudogap in the normal state. It 
is responsible for a lowering of Tc relative to the mean field 
value T*. CoUecting the quadratic terms in the above equa- 
tion, it follows that 

Tc = T*-— A2(Tc) . (47) 

Moreover, self consistency imposes a constraint on the mag- 
nitude of A(Tc) via 

A2(T) = f i^vl/e-'5^'[*lvI/2 / f i^vl/e-'^^I*! . (48) 
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FIG. 19 Pseudogap in the density of states above Tc in conventional 
superconductors, from lCohen et a/.(ll969l v\ = N{E) is the nor- 
malized density of states. 

It should be noted that Eq. ( 147 1 is consistent with the state- 
ment that A{Tc) and Tc lie on the BCS curve, since for small 
separation between Tc and T*, this curve is, below Tc, defined 
by 

AlcsiT)^^iT*-T). (49) 

The primary effect of fluctuations at the Hartree level is that 
pairing takes place in the presence of a finite excitation gap, 
reflecting the presence of very short lived, but pre-formed 
pairs. It should be clear that A is not to be confused with the 
superconducting order parameter Age, as should be evident 
from Fig. |3] We refer the reader to a more detailed discussion 
which is presented in Appendixicl 

Figure ^] illustrates how this pseudogap appears experi- 
mentally, via^^lotofthenomalized densit^ofstates_(Cohen 
et ai, l96^)T^i^W{E)'inthe normal state. This figure de- 
rives from tunneling measurements on Al-based films. Note 
the depression of N{Ep) at low energies. This is among the 
first indications for a pseudogap reported in the literature. 

Alternatively, fluctuations have also been discussed at the 
Gaussian (GqGq) level in which there is no need for self con- 
sistency, in contrast to the above picture. Here because the 
calculations are so m uch more tractable there have been very 
detailed applications (iLarkin and Varlamovll200 l') of essen- 
tially all transport and thermodynamic measurements. This 
Gaussian analysis led to unexpected divergences in the para- 
conductivity from the so-called Maki-Thompson term which 
then provided a motivation to go beyond leading order theory. 
Patton ( Patton, 1 97 1 a b ) showed how this divergence could be 
removed within a more self consistent GGq scheme, equiva- 



lent to self consistent Hartree theor y. Others argued dSchmidl 
Il969t iTucker and HalperirA 1 19711) that Hartree-Fock (GG) 
was more appropriate, although in this weak coupling, nar- 
row fluctuation regime, the differences between these latter 
two schemes are only associated with factors of 2. Neverthe- 
less with this factor of 2, A{Tc) and Tc are no longer on the 
BCS curve. It should be noted that a consensus was never 
fully reached ( H assing and Wilkins^ .1973.) by the community 
on this point. 



1 . Dynamics of Pre-formed Pairs Above Tc. Time Dependent 
Ginzburg-Landau Theory 

We turn again to Fig. |5] which provides a useful starting 
point for addressing bosonic degrees of freedom and their im- 
plications for experiment. It is clear that transport in the un- 
usual normal state, associated with the pseudogap (PG) phase 
(as well as thermodynamics) contains contributions from both 
fermionic and bosonic excitations. The bosons are not in- 
finitely long lived; their lifetime is governed by their inter- 
actions with the fermionSj^JnsomejnstancesJ^LarM and Var- 
lamov, 2001.) . the bosonic contributions to transport become 
dominant or even singular at Tc- Under these circumstances 
one can ignore the fermionic contributions except insofar as 
they lead to a lifetime for the pairs. 

To address the transport properties of these more well es- 
tablished pre-formed pairs, one appeals to a standard way of 
characterizing the dynamics of bosons - a phenomenological 
scheme known as time dependent Ginsburg- Landau (TDGL) 
theory. This theory has some microscopic foundations in dia- 
grammatic T-matrix approaches in which one considers only 
the Aslamazov-Larkin terms (which are introduced in Ap- 
pendix|n}. At the Gaussian level both TDGL and the T-matrix 
approaches are tractable. At the Hartree level things rapidly 
become more complicated, and i t is far easier to approach the 
problem ("Tan and L evinl |2004^ by adopting a strictly phe- 
nomenological TDGL. 

The generic Hartree-TDGL equation of motion for classical 
charged bosons interacting with electromagnetic fields ((/>, A) 
is given by 

/ d \ (-iV - e*A(x,<))^ 

- Hparr{T)ij{K,t)+D{x,t), (50) 

Here i/)(x, t) is a classical field variable representing the 
bosons which have vanishing chemical potential impair 3.1 Tc. 
The function £'(x, t) introduces a noise variable into the equa- 
tion of motion. Generally 7 is complex. 

We will explore the implications of this equation in our 
more quantitative discussion of Section lVi.EI 
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II. QUANTITATIVE DETAILS OF CROSSOVER 

A. T = 0, BEC Limit without Feshbach Bosons 

We begin by reviewing T = crossover theory in the BEC 
limit. Our starting point is the ground state wavefunction 
^0 of Eq. Q, along with the self consistency conditions of 
Eq. (|2}- For positive chemical potential ji, the quantity Asc(O) 
also corresponds to the energy gap for fermionic excitations. 
In the ground state the two energy scales A(0) and A^cCO) are 
degenerate, just as they are (at all temperatures) in strict BCS 
theory. 

It is convenient to rewrite these equations ^ in terms of the 
inter-fermion scattering length as 



E 

k 

E 



1 1 

2ek 2£'k 



1 



Sk 



T = 0. 



(51) 
(52) 



In the fermionic regime (/i > 0) these equations are essen- 
tially equivalent to those of BCS theory, although at weak 
coupling appropriate to strict BCS, little attention is paid to 
the number equation since /i = Ep is always satisfied. The 
more interesting regime corresponds to /i < where these 
equations take on a new interpretation. Deep inside the BEC 
regime it can be seen that 



AL(0) 



47r^2m|/x| 



(53) 



which, in conjunction with Eq. ( 1511 (expanded in powers of 

AL(0)^'): 



4m a. 



(2m) 



3/2 



1 



1 AL(0) 
16 



(54) 



yields 



To see the equivalence we associate the number of bosons 
npairs — n/2, the boson mass niB — 2m and the bosonic 
scattering length as — 2as. Here the bosonic chemical po- 
tential hb = 2/i + eg and we define eo — l/(2mag). This 
factor of 2 in as /as has received a fair amount of attention 
in the literature. Although it wiU be discussed later in the 
Review, it should be noted now that this particular numerical 
value is particular to the one channel description! Staiic et a/J 
2005 b), and differe nt from what is observed experimentally 
dPetrovef a/.Ll2004 . 

Despite these similarities with GP theory, the fundamental 
parameters are the fermionic Asc(O) and chemical potential 
/i. It can be shown that in this deep BEC regime the number 
of pairs is directly proportional to the superconducting order 
parameter 



ZoAlM 



where 



■m^as 



(57) 



(58) 



One may note from Eq. j57l i that the "gap equation" now cor- 
responds to a number equation (for bosons). Similarly the 
number equation, or constraint on the fermionic chemical po- 
tential defines the excitation gap for fermions, once the chem- 
ical potential is negative. In this way the roles of the two 
constraints are inverted ttRanderiaLll995h relative to the BCS 
regime. 

That a Gross-Pitaevskii approach captures the leading or- 
der physics at T = can also be simply seen by rewriting the 
ground_^ta te_WMefunction^_as_Qointe^^ bv Randeria ( Ran- 
deria, 1 19951) . Define v\./u\. = nir 



5*0 = const X n/j(l + 



const X exp 



|0) 



(59) 



2mai 



(55) Projecting onto a state with fixed particle number N yields 
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This last equation is equivalent (Perali et ai. 2003", Viverit 
et ai, 2004) to its counterpart in Gross-Pitaevskii theory. This 
theory describes true bosons, and is associated with the well 
known equation of state 



Anas 



(56) 



^ To make the equations simpler, in much of what we present in this re- 
view we will consider a contact potential interaction. The momentum 
sums are assumed to run to infinity. More generally, one can include a 
cutoff function (/j^ for the momentum integrals. The gap function thus 
becomes = A(^k. For the cuprate d-wave cuprate superconductors, 
Vk = cos k^a — cos kya. For our numerical work on atomic gases we 



use a Gaussian cutoff (/?^ 
numeiically feasible. 



'=^/'=o with fco taken to be 



as large as is 



const X ^y^??kc[.cl^.^ 



N/2 



|0) 



(60) 



This is a GP wavefunction of composite bosons with the im- 
portant proviso: that the characteristic size associated with 
the internal wavefunction 77k is smaller than the inter-particle 
spacing. It should be stressed that, except in the BEC asymp- 
totic limit, (where the "bosons" are non-interacting) there 
are essential differences between Eq. ( I60t and the Gross- 
Pitaevskii wavefunction. Away from this extreme limit, pair- 
pair interactions (arising indirectly via the fermions and the 
Pauli principle), as well as the finite size of the pairs, will 
differentiate this system from that of true interacting bosons. 
This same issue arises more concretely in the context of the 
collective mode spectrum of ultracold gases, as discussed in 
SectionlVCl 
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B. Extending conventional Crossover Ground State to 
T / 0: BEC limit without Feshbach Bosons 



How do we extend ( Chen et «/.L 11998*) this picture to fi- 
nite T? In the BEC limit, as shown in Fig. |3 fermion pairs 
are well established or "pre-formed" within the entire range 
of superconducting temperatures. The fundamental constraint 
associated with the BEC regime is that: for all T < Tc, there 
should, thus, be no temperature dependence in fermionic en- 
ergy scales. In this way Eqs. \51\ and ( 152 1 must be imposed 
at all temperatures T. 
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T <Tr 



(61) 
(62) 



It follows that the number of pairs at T = should be equal 
to the number of pairs at T = Tc. However, all pairs are con- 
densed at T = 0. Clearly, the character of these pairs changes 
so that at Tc, all pairs are noncondensed. To implement these 
physical constraints (and to anticipate the results of a more 
microscopic theory) we write 



^ condensed 



(T) 



noncondense\ 



(63) 
'^{T) (64) 



so that we may decompose the excitation gap into two contri- 
butions 



(65) 



where /S.sc{T) corresponds to condensed and Apg(T) to the 
noncondensed gap component. Each of these are proportional 
to the respective number of condensed and noncondensed 
pairs with proportionality constant Zq. At Tc, 



^ noncondensed 



pairs 



(66) 



where b{x) is the usual Bose-Einstein function and Vtq is the 
unknown dispersion of the noncondensed pairs. Thus 



o-^E 



b{^,,Tc) = -Z-\ (67) 



We may deduce directly from Eq. (I67> that A^^ — 
— t{Q), if we presume that below Tc, the noncondensed 
pairs have propagator 



t{Q) 



(68) 



It is important to stress that the dispersion of the pairs ilq 
cannot be put in by hand. It is not known a priori. Rather, 
it has be to derived according to the constraints imposed by 
Eqs. ( 16 U and ( 162^ We can only arrive at an evaluation of 
ilq after establishing the nature of the appropriate T-matrix 



theory. In what follows we derive this dispersion using the 
two channel model. We state the one channel results here 
for completeness. In the absence of FB, flq = B^q^, where 
Bo = l/(2AfQ). As a consequence of a "bosonic" self en- 
ergy, which arises from fermion-boson (as distinguished from 
boson-boson) interactions, the bosonic dispersion is found to 
be quadratic. This quadratic dependence persists for the two 
channel case as well. 



C. Extending Conventional Crossover Ground State to 
T / 0: r-matrix Scheme in the Presence of Feshbach 
Bosons 

To arrive at the pair dispersion for the noncondensed pairs, 
rig, we need to formulate a generalized T-matrix based 
scheme which is consistent with the ground state conditions, 
and with the T dependence of strict BCS theory. It is useful 
from a pedagog ical point to now include the effects of Fes- 
hbach bosons ( S taiic et al. 1 12004"). Our intuition concerning 
how true bosons condense is much better than our intuition 
concerning the condensation of fermion pairs, except in the 
very limited BCS regime. We may assume that A and /i 
evolve with temperature in such a way as to be compatible 
with both the temperature dependences of BCS and with the 
above discussion for the BEC limit. We thus take 



A(r) = -UeffY.^iT)'-^p^ 



E 



l-^ + 2^/(i?.) 



Sk 



(69) 



, (70) 



where Ey, = ^(ck - m)^ + A2(T). Note that in the pres- 
ence of Feshbach bosons, we use n to denote the density of 
fermionic atoms which is to be distinguished from the con- 
stituents of the Feshbach bosons. The total density of particles 
will be denoted by Equations ( I69t and 11701 will play a 
central role in this review. They have been frequently invoked 
in the literature, albeit under the presumption that there is no 
distinction between A(T) and Asc(T). 
Alternatively one can rewrite Eq. ( I69> as 
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(71) 



Clearly Eqs. \69\ and MQl are consistent with Eqs. ( I61> and 
(I62> since Fermi functions are effectively zero in the BEC 
limit. Our task is to find a T-matrix formalism which is com- 
patible with Eqs. ( I69> and ( I70> , and to do this we focus first on 
the noncondensed molecular bosons. Their propagator may be 
written as 



D{Q) 



1 



(72) 



We presume that the self energy of these molecules can be 
written in the form 



^B{Q) = -9\{Q)/[i + Uxm 



(73) 
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where xiQ) is as yet unspecified. This RPA-like self energy 
arises from interactions between the molecular bosons and the 
fermion pairs, and this particular form is required for self- 
consistency. 

noncondensed bosons in equilibrium with a condensate 
must necessarily have zero chemical potential. 



l^bo 



T <Tr 



(74) 



t = u + 



+ g 



Do 



Go 
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This is equivalent to the Hugenholtz-Pines condition that 

D-\0) = 0, T < Tc. (75) 
Using Eqs. <73> and j75> it can be seen that 



U-/j{0) + x{0) = 0, 



T <Tr 



(76) 



This equation can be made compatible with our fundamental 
constraint in Eq. \69\ provided we take 



x(Q) = ^G(i^)Go(g-i^), 



K 



where G{K) 



G{K)^[G^\K)-^{K)]- 



includes a self energy given by the BCS-form 

^{K) = -Ga{-K)^\ 



(77) 



(78) 



(79) 



which now involves the quantity A to be distinguished from 
the order parameter A^c- Note that with this form for Yi{K), 
the number of fermions is indeed given by Eq. ( 17m . 

More generally, we may write the constraint on the total 
number of particles as follows. The number of noncondensed 
molecular bosons is given directly by 



n,{T) = -Y,D{Q). 
Q 



(80) 



For T < Tc, the number of fermions is given via Eq. ( I70> . 
The total number (n*°*) of particles is then 



2ni, 



2nl 



n 



tot 



(81) 



where n° — 



is the number of molecular bosons in the 
condensate, and 0m is the order parameter for the molecular 
condensate. 

Thus far, we have shown that the condition that noncon- 
densed molecular bosons have zero chemical potential, can 
be made consistent with Eqs. ( I69> and (17 0> provided we con- 
strain x{Q) and S(-ft') as above. We now want to examine 
the counterpart condition on the fermions and their conden- 
sate contribution. The analysis leading up to this point should 
make it clear that we have both condensed and noncondensed 
fermion pairs, just as we have both condensed and noncon- 
densed molecular bosons. Moreover, these fermion pairs and 
Feshbach bosons are strongly admixed, except in the very ex- 
treme BCS and BEC limits. Because of noncondensed pairs. 



FIG. 20 Diagrammatic scheme for present T-matrix theory. Here U 
is the open-channel interaction, g is the coupling between open and 
closed channels. Go, G and Do are bare fermion, full fermion, and 
bare Feshbach boson Green's functions, respectively. U and goDo 
yield an effective pairing interaction, Ueff- 



we will see that the excitation gap is distinct from the super- 
conducting order parameter, although in the literature this dis- 
tinction has not been widely recognized. 

Just as the noncondensed molecular bosons have zero 
chemical potential below Tc we have the same constraint on 
the noncondensed fermion pairs which are in chemical equi- 
librium with the noncondensed bosons 



0, 



T <Tr 



(82) 



The T-matrix scheme in the presence of the Feshbach bosons 
is shown in Fig.|20l Quite generally, the T-matrix consists of 
two contributions: from the condensed (sc) and noncondensed 
or "pseudogap"-associated (pg) pairs. 



t = 



tsciQ) 



UeffiQ) 



l + Ucff{Q)x{Q)' 
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(83) 

^0, (84) 
(85) 



where we write A^c = A^c — g/pm, with Age = 
— [/^jj.(a_kiakT) and = (6q=o)- Here, the order pa- 
rameter is a linear combination of both paired fermions and 
condensed molecules. Similarly, we have two contributions 
for the fermion self energy 



^pg{K) ^Y.^{Q)Go{Q - K) , (86) 



where, as in BCS theory, 

^sc{K) ^Y.^srXQ)Go{Q - K) . 



(87) 



Without loss of generality, we choose order parameters A^^ 
and (/),„ to be real and positive with q < 0. Import antly, these 
two components are connected ( Kokkelmans et aU .200Z) by 
the relation = gAsc/[{i^ — 2/z)J7]. 

The vanishing of the pair chemical potential implies that 



t^„'{o) = u-)io) + x{o)^o, 



T<Tr_ 



(88) 



This same equation was derived from consideration of the 
bosonic chemical potential. Most importantly, we argued 
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above that this was consistent with Eq. ( I69> provided the 
fermion self energy assumes the BCS form. We now verify 
the assumption in Eq. (I79> . A vanishing chemical potential 
means that tpg{Q) is strongly peaked around Q = Q. Thus, 
we may approximate cMalv et al 1999a) Eq. ( I86> to yield 



where 



A2(T)^AL(r) + A2^(T). 



and we define the pseudogap A- 



pg 



(89) 



(90) 



(91) 



This bears a close resemblance to Eq. (18 0> . if we write the 
number of noncondensed fermion pairs and molecular bosons 

as 



^p{T) = ZAlJT) 



(92) 



Note that in the normal state (where iipair is nonzero) one 
cannot make the approximation of Eq. j89> . Referring back to 
our discussion of Hartree-approximated TDGL, a strong anal- 
ogy between Eq. ( I91> and ( I48> should be observed. There is 
a similar analogy between Eq. ( I69> and ( I49> : more details are 
provided in Appendix|C] We thus have that Eqs. ( 15 U and ( I89> 
with Eq. (|28} are alternative ways of writing Eqs. ( I69> and 
(I70> . Along with Eq. ( 19 U we now have a closed set of equa- 
tions for addressing the ordered phase. Moreover the propa- 
gator for noncondensed pairs can now be quantified, using the 
self consistently determined pair susceptibility. At moderately 
strong coupling and at small four- vector Q, we may expand to 
obtain 



tpg{Q) 



z- 



rig- -\- jj'pair ~t~ ^r^^ 

Consequently, one can rewrite Eq. < I9U as 

A2^(T) = Z-i^6(f^„T). 



(93) 



(94) 



Here we introduce a simple notation for the small Q expan- 
sions: 



with 



x(Q)-x(o) 
u;/f{Q)-u-/j{0) 



Zo 



Zo{i^r, 
Zg{i^r, 



2A2 



2^/(ek-M) 
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We have that Z = Zq + Zg, Zb = {Zq + Zg)/Zo, and the 
coefficient B = 1/(2M*) in fi, = Bq'^ is such that 



B 



BqZq + -^Zg 



The strong hybridization between the molecular bosons and 
fermion pairs should be stressed. Indeed, there is just one 
branch Q,q for bosonic-like excitations, so that at small Q 



D{Q) 



il — fig -f fJbboson + i^Q 



(95) 



where below Tc, ^ipair = fJ'boson = 0. Here Tq contains 
all the imaginary part of the inverse T-matrix and it vanishes 
rapidly as Q ^ 0. As the coupling is decreased from the 
BEC limit towards BCS, the character of the pairs changes 
from pre-dominantly molecular-boson like to (noncondensed) 
Cooper pairs, composed only of the fermionic atoms. 



D. Nature of the Pair Dispersion: Size and Lifetime of 
noncondensed Pairs Below Tc 



We may rewrite the pair susceptibility dChen et a/.Lll998i 
!1999^ of Eq. nil (after performing the Matsubara sum and 
analytically continuing to the real axis) in a relatively simple 
form as 



xiQ) 



E 

k 



/(i?k) - /(a-q) 

£^k-a-q + ^^ + ^o+' 



(96) 



where and are given by their usual BCS expressions in 
terms of A and ^ = Ck ^ In the long wavelength, low 
frequency limit, the inverse of tpg can be written as 



aifi^ + Zo(fi 



q- 



2AP 



+ fJ-pair + irg). 



(97) 



We are interested in the moderate and strong coupling 
cases, where we can drop the aiil^ term in Eq. i91i . and hence 
we have Eq. ( I93l l with 



fir 



2M* 



(98) 



This establishes a quadratic dispersion and defines the effec- 
tive pair mass, M*. Analytical expressions for this mass are 
possible via a small q expansion of x, in Eq- ( I96> . It is impor- 
tant to note that the pair mass reflects the effective size £,pg of 
noncondensed pairs. This serves to emphasize the fact that the 
q^ dispersion derives from the compositeness of the "bosons", 
in the sense of their finite spatial extent. A description of the 
system away from the BEC limit must accommodate the fact 
that the pairs have an underlying fermionic character. This 
pair mass has a different origin from the mass renormaliza- 
tion associated with real interacting bosons. There one finds a 
mass shift which comes from a Hartree approximate treatment 
of their g-dependent interactions. Finally, we note that £^pg is 
comparable to the size ^ of pairs in the condensate. In the 
weak coupling BCS limit, a small Q expansion of xo shows 
that, because the leading order term in 57 is purely imaginary, 
the contribution cannot be neglected. The T-matrix, then, 
does not have the propagating q^ dispersion, discussed above. 
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FIG. 21 Typical behavior of Tc in the homogeneous case. For the 
purposes of illustration, this calculation includes Feshbach effects. 
Tc follows BCS for large magnetic detuning parameter and ap- 
proaches the BEC asymptote 0.218 for large negative uq. It reaches a 
maximum at unitarity, and has a minimum when the fermionic chem- 
ical potential changes sign. 



It follows from Eq. ( I96> that the pair lifetime 



^0 

+ [/(i^k) - /(Ck-q)] V15{E^ - Ck-q + n) . 



(99) 



Here Tq reflects the rate of decay of noncondensed bosons 
into a bare and dressed fermion. Note that the excitation gap 
in significantly restricts the contribution from the 6 func- 
tions. Thus, the decay rate of pair excitations is greatly sup- 
pressed, due to the excitation gap (for fermions) in the super- 
conducting phase. Even in the normal phase pairs live longer 
than one might have anticipated from "Pauli blocking" argu- 
ments which are based on evaluating Tq in a Fermi liquid 
state. However, the same equations are not strictly valid above 
Tc, because Eq. ( I89> no longer holds. To compute the pair life- 
time in the normal state requires a more extensive calculation 
involving the full T -matrix self consistent equations and the ir 
numerical solution jJanko et a/.Lll997HMalv ef aail999albh . 



E. Tc Calculations: Analytics and Numerics 

Calculations of Tc can be performed using Eqs. \69\ and 
j70l l along with Eq. j91l l. The T-matrix is written in the ex- 
panded form of Eq. ( I93> . which is based on the pair dispersion 
as derived in the previous section. For the most part the cal- 
culations proceed numerically. 

A typical curve is plotted in Fig. |^ where the three 
regimes BCS, BEC and PG are indicated. For positive, but de- 
creasing the magnetic detuning vq, Tc follows the BCS curve 
until the "pseudogap" or A(Tc) becomes appreciable. After 
its maximum (near the unitary limit, at slightly negative vq), 
Tc decreases, (as does /i), to reach a minimum at /i w 0. This 
decrease in Tc reflects the decreasing number of low energy 
fermions due to the opening of a pseudogap. Beyond this 
point, towards negative vq, the system is effectively bosonic, 
and the superconductivity is no longer hampered by pseudo- 
gap effects. In the presence of FB, the condensate consists 



of two contributions, although the weight of the fermion pair 
component rapidly disappears. Similarly Tc rises, although 
slowly, towards the ideal BEC asymptote, following the in- 
verse effective boson mas s. The co rresponding curve based 
on the NSR approach (Nozieres and Schmitt-Rink. 1985) has 
only one extremum , but nevertheless the overall magnitudes 
are not so different ( iMilstein et a /.L 120021: lOhashi and GriffiiiL 
|2002). 

Analytic results are obtainable in the near-BEC limit only. 
The general expression for the effective mass 1/A/q in this 
limit is given by 



M* 



(100) 



where here FB effects have been dropped for simplicity. In 
what follows, we expand Eq. (I100> in powers of na^ and ob- 
tain after some algebra 



Mn* w 2m 1 + 



(101) 



We now invoke an important constraint, derived earlier in 
Eq. (I66> which corresponds to the fact that at Tc all fermions 
are constituents of uncondensed pairs 

=Y,h{n^,Tc) . (102) 



From the above equation it follows that (M*Tc)'^/^ 
const, which, in conjunction with Eq. jlOU implies 

Tr - TP 7ra?n 



cx n 



TP 



(103) 



Here T^ is the transition temperature of the ideal Bose gas 
with Mq = 2m. This downward shift of Tc follows the ef- 
fective mass renormalization, much as expected in a Hartree 
treatment of GP theory at Tc. Here, however, in contrast to 
GP theory for a homogeneous system with a contact potential 
llDalfovo et a/.L [1999). there is a non-vanishing renormaliza- 
tion of the effective mass. This is a key point which under- 
lines the importance in this approach of the fermionic degrees 
of freedom, even at very strong coupling. 



III. SELF CONSISTENCY TESTS 

All T-matrix approaches to the BCS-BEC crossover prob- 
lem must be subject to self consistency tests. Among these, 
one should demonstrate that normal state self energy effects 
within a T-matrix scheme are not associated w ith superflu- 
idity or superconductivity ("Ko sztin et a/.LfT99^ . While this 
seems at first sight straightforward, all theories should be put 
to this test. Thus for a charged system, a meaningful result 
for the superfluid density ps requires that there be an exact 
cancellation between diamagnetic and paramagnetic current 
contributions at Tc. In this way self energy effects in the num- 
ber equation and gap equation must be treated in a consistent 
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fashion. In the following we address some of these self con- 
sistency issues in the context of the asymmetric (x ~ GGq) 
BCS-BEC crossover scheme. Other crossover theories with 
different forms for the pair susceptibility should also be sub- 
iected to these tests. We begin with a surn mary of calculations 
(iKosztin et oil 120001: Istaiic et adl2003ah of the Meissner ef- 
fect, followed by a treatment of the collective mode spectrum. 
Our discussion in this section is presented for the one channel 
problem, in the absence of Feshbach effects. 



A. Important check: behavior of ps 

The superfluid density may be expressed in terms of the 
local (static) electromagnetic response kernel A'(0) 



r^(0) 



(104) 



with the current-current correlation function given by 

Pa/3(Q) = -2e2^A„(i^,i^ + Q)G(/^ + Q) 



K 



(105) 



Here the bare vertex X{K, K 
consider Q — (q, 0), with q 



Q) = ;i^(k + q/2) and we 
> 0. We write A^*^ = A + 



(5Apg + 5A.SC, where the pseudogap contribution (5Apg to the 
vertex correction will be shown in Appendix IbI to satisfy a 
Ward identity below Tc 



dk 



(106) 



By contrast for the superconducting contributions, one has 



dk 



(107) 



This important difference in sign is responsible for the fact 
that the Meissner effect is associated with superconductivity, 
and not with a normal state self energy. 

The particle density n, after partial integration can be 
rewritten as n = —(2/3) J2k^ ' 9kG{K). Then, as a re- 
sult of Dyson's equation, one arrives at the following general 
expression which relates to the diamagnetic contribution 



m 



dk 



k- 



dk 



(108) 

Now, inserting Eqs. ( I108> and ( I105> into Eq. ( I104> one can 

see that the pseudogap contribution to iis drops out by virtue 
of Eq. ( I106t : we find 



= ^J2^HK)k-(sA 



K 



dk 



(109) 



We emphasize that the cancellation of this pseudogap con- 
tribution to the Meissner effect is the central physics of this 



analysis, and it depends on treating self energy effects in a 
Ward-identity- consistent fashion. 
We also have that 

6AsciK+Q,K) = A''^^Goi-K-Q)Goi-K)X{K+Q,K) . 

(110) 

Inserting Eqs. (ETJl, and ( fTTOt into Eq. ( fT09l . after cal- 
culating the Matsubara sum, one arrives at 



l-2/(^k) 



2Eu 



(111) 



This expression can be simply rewritten in terms of the BCS 
result for the superfluid density 



V m / V TO / 



BCS 



(112) 



Here [ns/m)-^'-'^ is just (ns/m) with the overall prefactor 
A^^ replaced with A^. 

Finally we can rewrite Eq. (II 12> using Eq. ( I90> as 



(-) 



A2 

pg 

A2 



(-) 

V TO / 



BCS 



(113) 



In this form it is evident that (via A^^) pair excitations out of 
the condensate are responsible for a suppression of the super- 
fluid densi ty relative to that o btained from fermionic excita- 
tions only \Staiic et alJ \2003l\) . 



B. Collective Modes and Gauge Invariance 

The presence of pseudogap self energy effects greatl y com- 
plicat es the computation of collective modes (Kosztin et al\. 
I2000h . This is particularly apparent at nonzero temperature. 
Once dressed Green's functions G enter into the calculational 
schemes, the collective mode polarizabilities and the elec- 
tromagnetic (EM) response tensor must necessarily include 
vertex corrections dictated by the form of the self-energy S, 
which depends on the T-matrix which, in turn depends on the 
form of the pair susceptibility x- These necessary vertex cor- 
rections are associated with gauge invariance in the same way, 
as was seen for ps, and discussed in AppendixlHl Collective 
modes are important in their own right, particularly in neutral 
superfluids, where they can be directly detected as signatures 
of long range order They also must be invoked to arrive at 
a gauge invariant formulation of electrodynamics. It is rela- 
tively straightforward to introduce these collective mode ef- 
fects into the electromagnetic response in a completely gen- 
eral fashion that is required by gauge invariance. The diffi- 
culty is in the implementation. 

In the presence of a weak externally applied EM field, with 
four-vector potential — {(j>, A), the four-current density 
= {p, J) is given by 



jf^{Q) = i^^-(Q)A,(g) , 

where K^^'^{Q) is the EM response kernel. 



(114) 
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The incorporation of gauge invariance into a general mi- 
croscopic theory may be implemented in several wa ys. Here 
we do so via a generalized matrix Kubo formula ( Kul ik et all 
[198 1 ) in which the perturbation of the condensate is included 
as additional contributions Ai + i/S.2 to the applied external 
field. These contributions are self consistently obtained (by 
using the gap equation) and then eliminated from the final ex- 
pression for K^^'^ . We now implement this procedure. Let 
ryi,2 denote the change in the expectation value of the pair- 
ing field f}i 2 corresponding to Ai 2. For the case of an s- 
wave pairing interaction U < Q, the self-consistency condi- 
tion A1.2 = ?/7?i.2/2 leads to the following equations: 



2Ai 



m 



2A2 

\u\ 



R^^A, + OiiAi + Q12A2 , (115b) 
R^^'A, + Q21A1 + Q22A2 , (115c) 



where 



ifo^''(a;,q) - P^^{uj, q) + — - 5^°) (116) 

m 

is the usual Kubo expression for the electromagnetic re- 
sponse. We define the current-current correlation function 

P'^"^(T,q) = -*e(T)([j^(r,q),j-(0,-q)]). In the above 
equation, g^'' is a (diagonal) metric tensor with elements 
(1,-1,-1,-1). We define 

i?^^(r,q) = -*0(r)([j'^(r,q),,7,(O,-q)]) (117) 
with /i = 0, . . . , 3, and i, j = 1,2; and 

Q„(r,q) - -je(T)([77,(T,q),^,(0,-q)]) (118) 
Finally, it is convenient to define 

Q^^=2/U + Q^^. (119) 

In order to demonstrate gauge invariance and reduce the 
number of component polarizabilities, we first rewrite iiT^" in 
a way which incorporates the effects of the amplitude contri- 
butions via a renormalization of the relevant generalized po- 
larizabilities. 



^0 — -"-0 



Qii 



(120) 



It can be shown, after some analysis, that the gauge invari- 
ant form for the response tensor is given by 



K"" = K 



n ,K'^ n 



(121) 



The above equation satisfies two important requirements: it 
is manifestly gauge invariant and, moreover, it has been re- 
duced to a form that depends principally on the four-current- 
current correlation functions. (The word "principally" appears 



because in the absence of particle-hole symmetry, there are 
effects associated with the order parameter amplitude contri- 
butions that add to the complexity of the calculations). 

The EM response kernel of a superconductor contains a 
pole structure that is related to the underlying Goldstone bo- 
son of the system. Unlike the phase mode component of the 
collective mode spectrum, this Anderson-Bogoliubov (AB) 
mode is independent of Coulomb effects. The dispersion of 
this amplitude renormalized AB mode is given by 



(122) 



For an isotropic system K'q — K'q Sap, and Eq. J122> can 
be rewritten as 



2 7,w00 







c^^K'f - 2..q„if'°" = , 



(123) 



with a = 1, 2, 3, and in the last term on the left hand side 
of Eq. MTil a summation over repeated Greek indices is as- 
sumed. 

It might seem surprising that from an analysis which incor- 
porates a complicated matrix linearized response approach, 
the dispersion of the AB mode ultimately involves only 
the amplitude renormalized four-current correlation func- 
tions, namely the density-density, current-current and density- 
current correlation functions. The simplicity of this result is, 
nevertheless, a consequence of gauge invariance. 

At zero temperature K'^" vanishes, and the sound-like AB 
mode has the usual linear dispersion lo = ujq_ = c|q| with the 
"sound velocity" given by 



(124) 



We may, thus, interpret the AB mode as a special type of col- 
lective mode which is associated with A^, = 0. This mode cor- 
responds to free oscillations of Ai 2 with a dispersion oj — cq 
given by the solution to the equation 



detlQy I = Q11Q22 - Q12Q21 = 



(125) 



The other collective modes of this system are derived by in- 
cluding the coupling to electromagnetic fields. Within self- 
consistent linear response theory the field 6(p must be treated 
on an equal footing with Ai 2 and formally can be incorpo- 
rated into the linear response of the system by adding an extra 
to the right hand side of Eqs. 



TT5ab . mm . and 
il 15c> . Note that, quite generally, the effect of the "external 
field" S(p amounts to replacing the scalar potential A'^ — cf) hy 
A'^ ^ <f> = </) + Scj). In this way one arrives at the following set 



term Klf' 



of three linear, homogeneous equations for the unknowns 1 
Ai, and A2 



= i?i"<50 + QiiAi 
= R^'^5<p + Q2iAi 



Q12A2 , (126a) 
Q22A2 , (126b) 



i?°iAi +i?''2^2 , (126c) 



where V is an effective particle-hole interaction which may 
derive from the pairing channel or, in a charged superconduc- 
tor, from the Coulomb interaction. The dispersion of the col- 
lective modes of the system is given by the condition that the 
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FIG. 22 Temperature dependence of the real and imaginary parts of 
tfie AB collective mode velocity for moderate (solid lines) and weak 
(dashed lines) coupling. The velocity vanishes at Tc for both cases. 



above equations have a nontrivial solution 



l/U 



Ql2 

Q22 + i/U 



^0 



1/V 



(127) 



In the BCS limit where there is particle-hole symmetry Q12 = 
Q21 — — R'^^ — and, the amplitude mode decouples 
from the phase and density modes; the latter two are, however, 
in general coupled. 

The above formalism (or its equivalent RPA variations 
(Belkhir and Randeria, 1992; Cote and Griffin, 1993)) has 
been applied to address collective modes in the crossover sce- 
nario. The most extensive studies have been at T = based 
on the ground state of Eq. Q. There one finds a smooth 
change in the character of the Anderson-Bogoliubov (AB) 
mode. At weak coupling one obtains the usual BCS value 
c = vp/^/3. By contr ast at strong coupUn^, the collective 
mode spectrum reflects ( Belkhir and Randeria, '1992') an effec- 
tive boson-boson interaction deriving from the Pauli statistics 
of the constitu ent fermions. This is most clearly seen in jel- 
lium models ( Griffin and Ohashil Eool iKosztin ef"aZl l200(i) 
where the AB sound velocity is equivalent to that predicted 
for a 3D interacting Bose gas c = [iTrnaB/m^]^/^. Here, 
as earlier, the inter-boson scattering length is twice that of the 
inter-fermionic counterpart, at strong coupling.^ 

The effects of finite tem perature on the AB mode have been 
studied in (.Kosztin et aU 2Q0Q) . by making the approxima- 
tion that the temperature dependence of the order parameter 
amplitude contribution is negligible. Then the calculation of 
this dispersion reduces to calculations of the electromagnetic 



response, as discussed, for example in pg. See also Appendix 
151 At finite T, the AB mode becomes damped and the real and 
imaginary parts of the sound velocity have to be calculated 
numerically. Here one finds, as expected, that the complex 
c ^ as T ^ Tc. The results are plotted in Fig.|22]for both 
moderate and weak coupling. Because the real and imaginary 
parts are comparable there, the mode is highly damped near 



C. Investigating the Applicability of a Nambu Matrix 
Green's Function Formulation 

Diagrammatic schemes appropriate to BCS superconduc- 
tors are based on a Nambu matrix Green's function approach. 
The off-diagonal or anomalous Green's functions in this ma- 
trix are given by 



F{K) = A,cG{K)Go{^K) 



(128) 



as well as its Hermitian conjugate. This Nambu formalism 
was designed to allow study of perturbations to the BCS 
state, due to, for example, external fields or impurities. For 
these perturbations the central assumption is that they act on 
both the "normal" (with BCS Green's function G{K) from 
Eq. (I30» and anomalous channels in a symmetrical way. 

Understanding as we now do that BCS theory is a very 
special case of superconductivity, this raises the caution that 
once one goes beyond BCS, some care should be taken to jus- 
tify this Nambu approach. At the very least the distinction 
between the order parameter A^c and the excitation gap A 
raises ambiguity about applying a diagrammatic theory based 
on Nambu-Gor'kov Green's functions. 

More importantly, in treating the effects of the pseudogap 
self energy Epg as we do here, it should be clear that this self 
energy does not play a symmetric role in the anomalous and 
normal channels. It is viewed here as an entirely normal state 
effect. However thi s theory accommodates the equivalent 
(iKosztin et al\. 12000) of the Nambu-Gor'kov "F"- function 
in general response functions such as in the Maki-Thompson 
diagrams of Appendix 151 through the asymmetric combina- 
tion GGo which always arises in pairs (for example, as the 
FF combinations of BCS theory). In this regard the GGq 
formalism appears to differ from all other T-matrix schemes 
which are designed to go below Tc, in that the Nambu scheme 
is not assumed at the start. Nevertheless, many features of 
this formalism seem to naturally arise in large part because 
of Eq. ( I128> . which demonstrates an intimate connection be- 
tween GGo and the conventional diagrams of BCS theory. 



IV. OTHER THEORETICAL APPROACHES TO THE 
CROSSOVER PROBLEM 



* In the neutral case, for the full collective modes of Eq. 11271 . there are nu- 
merical differences (of order unity) in the prefactors of the mode frequency, 
so that the collective modes of the crossover theory are not strictly the same 
as in GP theory. In these calculations V in Eq. 11271 should be associ- 
ated with the pairing interaction, which enters in the context of fermionic 
density-density correlation effects within the particle-hole channel. 



A. T / Theories 

We have noted in Section lLGl that there are three main theo- 
retical approaches to the crossover problem based on T-matrix 
theories. Their differences are associated with different forms 
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for the pair susceptibility x- The resulting calculations of 
show similar variations. When two full Green's functions are 
prese nt in y (as in the FL EX approach), Tc varies monotoni- 
cally (Haussmann, 1993) with increasing attractive coupling, 
approaching the ideal gas Bose-Einstein asymptote from be- 
low. When two bare Green's functions are present in %, as 
in the work of Nozieres and Schmitt-Rink, and of Rande- 
ria and c o-workers, then Tc o vershoots (Sa de M elo e t ai, 
[^93; No zieres and Schmitt-R ink. 1985; Randerial ll995l) the 
BEC asymptote and ultimately approaches it from above. Fi- 
nally when there is one bare and one dressed Green's function, 
Tc first overshoots and then decreases to a minimum around 
M = 0, eventua lly approaching the asymptote from below 
dMalv et al. I[m9b|). This last appears to be a combination 
of the other two approaches. Overall the magnitudes are rela- 
tively similar and the quantitative di fferences are small. There 
is, fin ally, another body of work ( iGiovannini and Berthodl 
1200 Ih based on the GGq scheme which should be noted here. 
The approach by the Swiss group seeks to make contact with 
the phase fluctuation scenario for high Tc superconductors and 
addresses STM and ARPES data. 

Bigger differences appear when these T-matrix theories are 
extended below Tc- Detailed studies are most extensive for 
the NSR approach. Below Tc one presumes that the T-matrix 
(or x) contains only bare Green's functions, but these func- 
tions are now taken to correspond to their Nambu matrix form, 
with the order parameter Age appearin g in the dispersi on rela - 
tion for excited fermions (Fi eri et all 12004: Randeria. ll995l) . 
A motivation for generalizing the below- Tc T-matrix in this 
way is that one wants to connect to the collective mode spec- 
trum of the superconductor, so that the dispersion relation for 
pair excitations is Vlq k, cq. In this way the system would be 
more dkectly analogous to a true Bose system. 

Self energy effects are also incorporated below Tc but only 
in the number equation, either in the approximate manner of 
NSR (Griffin and Ohashi, 2003; Rande ria, 19 95) o r thro ugh 
use of the full Dyson resummation dPieri et all l2004l) of 
the diagonal Nambu-matrix component G(Eo). For the lat- 
ter scheme, Strinati's group has addressed pseudogap effects 
in some detail with emphasis on the experimentally observed 
fermionic spectral function. Some concern can be raised that 
the fermionic excitations in the gap equation do not incor- 
porate this pseudogap, although these pairs are presumed to 
emerge out of a normal state which has a pseudogap. Indeed, 
this issue goes back to the original formulation of NSR, which 
includes self energy effects in the number equation, and not in 
the gap equation. 

Tchernyshyov dTchernvshvovl fl997h presented one of the 
first discussions below Tc for the FLEX-based T-matrix 
scheme. He also addressed pseudogap effects and found a 
suppression of the fermion density of states at low energy 
which allows for long-lived pair excitations inside this gap. 
At low momenta and frequencies, their dispersion is that of a 
Bogoliubov-sound-like mode with a nonzero mass. Some of 
these ideashave^eenrecentl^exte^ and Ran- 

nmger, 20Q3F to apply to the boson-fermion model. 

The work of Micnas and collaborators dMicnas efa/.lll995h 
which is also based on the FLEX scheme predates the work 



of Tchernyshyov, although their initial focus was above Tc in 
two dimensional Hubbard models. Here the fermionic spec- 
tral functions were studied in detail, using a numerically gen- 
erated solution of the T-matrix equations. A nice review of 
their large body of work on "real space pairing" in the cuprates 
was p resented based on a variety of diffe rent model Hamilto- 
nians ("Micnas and RobaszkiewiczVl998"). 

An extensive body of work on th e FLEX scheme has been 
contributed by Yamada and Yanatse dYanase et a/.Ll2003h both 
above and below Tc- They point out important distinctions 
between their approach and that of NSR. The effects of the 
broken symmetry are treated in a gene rahzed Nam bu formal- 
ism, much as assumed by Fieri et al. (('Fierief a/.','2004)) and 
by Griffin and Ohashi ((Griffin and Ohashi, 2003)), but here 
the calculations involve self energy effects in both the number 
and gap equations. Their work has emphasized the effects of 
the pseudogap on magnetic properties, but they have discussed 
a wide variety of experiments in high Tc and other exotic su- 
perconductors. 

The nature of the ground states which result from these two 
alternatives (x ~ GG and x ~ GqGo ) has yet to be clearly 
established. It should be noted that, despite the extensive body 
of work based on the NSR or GqGq scheme, no detailed pic- 
ture has been presented in the literature of the nature of the 
ground state which results when their calculations at Tc are 
extrapolated dPieri et al 1 12004 down to T = 0. In large 
part, the differences between other work in the literature and 
that summarized in Sec tions fll. Cl are due to the whether (as in 
NSR-based papers) or not (as here) the superconducting order 
parameter A^c alone characterizes the fermionic dispersion 
below Tc. 

One might well ask the question: because the underlying 
Hamiltonian [Eq. dl IH is associated with inter-fermionic, not 
inter-bosonic interactions, will this be reflected in the near- 
BEC limit of the crossover problem? The precise BEC limit 
is, of course, a non-interacting, or ideal Bose gas, but away 
from this limit, fermionic degrees of freedom would seem 
to be relevant in ways that may not be accounted for by the 
analogue treatment of the weakly interacting Bose gas. The 
most extensive study (albeit, above Tc) o f this issue is due to 
Fieri and Strinati (tFieri and Strinatiil2000<) . Their work, im- 
portantly, points out the inadequacies of T-matrix schemes, 
particularly at strong coupling. While the ground state in their 
calculations is unknown, it is necessarily different from the 
conventional crossover state of Eq. (^. There is much in- 
tuition to be gained by studying this simplest of all ground 
states, as outlined in this Review, but it will clearly be of great 
value in future to consider states, in which, for example, there 
is less than full condensation. 

Finally, we note that studies of non-s wave orbital pairing 
states in the context of the BCS-BEC crossover problem have 
rapidly proliferated in the literature. This is, in large part be- 
cause of the observed d-wave pairing in the cuprates. Indeed, 
Leggett's initial work (Leasett, 198(3), a ddressed p-wave pair- 
ing in "^He. Randeria and co-workers ('Randeria','l995) also 
considered non-s wave p airing. Studies (Andrenacci et ail 
1 19991 : Iden Herto^ Il999h of this crossover for a d-wave lat- 
tice case are plentiful at T = 0. Even earlier work was done 
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dChen et al. Ill999> on extensions of this ground state to in- 
clude nonzero temperature. Moreover, a d -wave c ontinuum 
model has also been considered (Ouintanill a et a/.L 12002). It 
is expected that in the near future cold atom experiments will 
begin to address I ^ pairing via appropriat e choices of the 
Feshbach resonances llHo and Zaharievil2004) . 



B. Boson-Fermion Models for High Superconductors 

The notion that it might be appropriate to include bosonic 
degrees of freedom in high Tc superconductors g oes back to 
T.D. Lee and R. Friedberg who introduced ( Friedber g and Leei 
Il989albh the "boson-fermion" model within three years of the 
discovery of cuprate superconductivity. The work of Lee and 
colleagues is based on the Hamiltonian of Eq. (II 1> without 
the direct fermion-fermion interaction U . To our knowledge 
this is the first study which argues for something akin to "real 
space pairing" in these materials. Interestingly, this work was 
motivated only by the observed short coherence length. Pseu- 
dogap effects, which make this case much stronger, had not 
yet been properly characterized. The analogies with He"* and 
a more BEC-like picture of superconductivity were discussed, 
as well as the superfluid density, and energy gap structure, al- 
beit for s-wave pairing. 

Subsequent work on this model Hamiltonian in the con- 
text of the high Tc copper oxides was pursued by Ranninger 
and Micnas and their collaborators. Indeed, Ranninger and 
Robaszkiewicz (Ranninger and Robaszkiewicz, 1985) were 
the first to write down the boson-fermion model, although 
initially, in an entirely different context, the bosons were as- 
sociated with phonons rather than electron pairs, and they 
were presumed to be localized. In this context they addressed 
mean field theoretic solutions be low Tc as well as more re- 
cently dPomanski and Ranningeii l2003) an approach based on 
renormalization group techniques. The physical picture which 
emerges leads to two gap structures below T^, one associated 
with a pseudogap and another with the superconductivity. In 
this way it appears to be very different from approaches based 
on Eq. ([D- 



V. PHYSICAL IMPLICATIONS: ULTRACOLD ATOM 
SUPERFLUIDITY 

A. Homogeneous case 

In this section we summarize the key characteristics of 
fermionic superfluidity in ultracold gases, in the homogeneous 
situation, without introducing the trap potential. Our results 
are based on numerical solution of the coupled Eqs. ( I69> . (I70> 
and (|9T)- The upper left panel of Fig.|23]plots the fermionic 
chemical potential, the Feshbach boson condensate ratio and 
the inverse scattering length as a function of v^. Here we have 
chosen what we believe is the physically appropriate value 
for the Feshbach coupling go = —^OEp/kp for "'"K, as de- 
termined from the experimentally measured scattering length 
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FIG. 23 (Color) Characteristics of the ground state and as 
a function of the magnetic detuning in a homogeneous case. 
Shown are the chemical potential /i, fraction of molecular conden- 
sate 2n(,o inverse scattering length m/47ras, gap A(0), inverse 
pair size m/A-K^, dimensionless relative pair size kp^ at T = as 
well as Tc and A{Tc). The pair size is roughly the same as inter- 
particle spacing, fcf^ ~ 1, at unitarity. As z/q decreases, the gap A 
increases and reaches a maximum, and then decreases slowly and ap- 
proaches its BEC asymptote (not shown) from above. Here we take 
go = —AOEp/k'^^ and Uo — —SEp/kp, and we use the conven- 
tion for units Ep — kp — h = ks = 2m = 1. A large exponential 
momentum cutoff, k^/kp — 80 is assumed in the calculations. 

data and resonance width. ^ This parameter is roughly an order 
of magnitude larger for ^Li. Here i'q, the gaps and /i are all in 
units of Ep, and the plots, unless indicated otherwise, are at 
zero temperature. 

The upper right panel plots the excitation gap at Tc as well 
as the gap at T = 0. The lower left hand panel indicates Tc 
along with the inverse pair size in the condensate. Finally 
the lower right panel plots ^ itself, along with . 

One can glean from the figure that the Feshbach boson frac- 
tion decreases, becoming negligible when the chemical po- 
tential passes through zero. This latter point marks the on- 
set of the PG regime, and in this regime the condensate con- 
sists almost entirely of fermionic pairs. The upper limit of the 
PG regime, that is, the boundary line with the BCS phase, is 
reached once the pseudogap, A(Tc) is essentially zero. This 
happens when /i is close to its saturation value at Ep. At this 
point the pair size rapidly increases. 

These results can be compared with those derived from 
a smaller value of the Feshbach coupling constant go — 
— lOEp/k^^ shown in Fig. 231 Now the resonance is effec- 



' We note that the precise values of go and Uo are not critical for either ®Li or 
^"K as long as they give roughly the right (large) resonance width AB and 
background scatter length ai,g. In fact, these parameters usually change 
during experiment since the total particle number A'^ (and hence Ep and 
kp) decreases in the process of evaporative cooling. In order to compare 
with experiment, one needs to compare at the coiTect value of kpas- 
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FIG. 24 (Color) Characteristics of the ground state and T^, as in 
Fig. l23l except for a smaller go. In comparison with Fig. 1231 the Fes- 
hbach resonance becomes much narrower. Here go — —lOEp /k^^ 
and Uo = ~3EF/k%. 



FIG. 26 (Color) Characteristics of the ground state and Tc, as in 
Fig- mi except now Uo = 0. As a consequence of the vanishing of 
Uo, the gap A is monotonic and approaches its BEC asymptote from 
below as uo decreases. Here go = — WEp/Kp^^. 




FIG. 25 (Color) Characteristics of the ground state and Tc in the 
one-channel model, as in Fig.|23|but as a function of the attraction 
strength U. Here go = 0. Different from the two-channel shown in 
Figs. l23l24l case. the gap A increases indefinitely with \ U\. The units 
for U are Ep/k^. Note that \U\ increases to the left. 



0.3 




FIG. 27 (Color) Behavior of Tc (main figure) and A(Tc) (inset) at 
fixed go = —^QEp/ky.^ with variable Uo- As Uo — oo, unitarity 
is approached from the BCS side at large positive vo- At the same 
time, A oc \g\ also decreases. 



lively narrower. Other qualitative features remain the same as 
in the previous figure. 

One can plot the analogous figures in the absence of Fes- 
hbach effects. Here the horizontal axis is the inter-fermionic 
interaction strength |[/|, as is shown in Fig. 1251 Three essential 
differences can be observed. Note, first the obvious absence of 
the Feshbach boson or molecular condensate ribo- There is, of 
course, a condensate associated with pairs of fermions (Asc) 
and these pairs will become bound into "fermionic molecules" 
for sufficiently strong attraction. Secondly, note that the exci- 
tation gap A(0) is monotonically increasing as \U\ increases 
towards the BEC limit. By contrast, from the upper right 
panels of Fig.|24]one can see that when Feshbach effects are 



present A(0) « A(Tc) decreases towards zero in the extreme 
BEC limit. It can also be seen that the shape of the scattering 
length curve vs U is different from the plots in the previous 
two figures. Here Os more rapidly increases in magnitude on 
either side of the unitary limit. Nevertheless it is important to 
stress that, except for the nature of the Bose condensate in the 
BEC limit, the physics of the Feshbach resonance-induced su- 
perfluidity is not so qualitatively different from that associated 
with a direct inter-fermionic attraction. 

Figure |26l presents the analogous plots for the case when 
the direct attraction between fermions U vanishes. Here the 
coupling constant is the same as in Fig. |23 and pseudo- 
gap effects are present here just as in all other cases. The 
most notable difference in the physics is that the formation 
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FIG. 28 Typical density of states (DOS) N{lu) vs energy uj in the 
three regimes at indicated temperatures. There is no depletion of the 
DOS at Tc in the BCS case. By contrast, in the BEC case, there is 
a large gap in DOS at Tc, roughly determined by In the inter- 
mediate, i.e., pseudogap, regime, there is a strong depletion of the 
DOS already at Tc. Once T decreases below Tc, the DOS decreases 
rapidly, signalling the superfluid phase transition. 



of the Cooper condensate is driven exclusively by coupling 
to the molecular bosons. A contrast between the two fig- 
ures is in the behavior of A{Tc) and A(0). For the U = 
case and in the BEC limit, A(0) approaches the asymptotic 
value i \g\ y^n^) from below. When the system is the fermionic 
regime (/^ > 0), aside from quantitative details, there is very 
little to distinguish the effects of U from those of g, except 
when these parameters are taken to very extreme limits. 

Finally, in Fig. |^we plot vs i^o in the presence of Fes- 
hbach effects with variable Uq, from weak to strong back- 
ground coupling. The lower inset shows the behavior of the 
excitation gap at Tc. With very strong direct fermion attrac- 
tion U, we see that Tc has a very different dependence on i/q. 
In this limit there is a molecular BEC to PG crossover (which 
may be inaccessible in actual experiments, since U is not suf- 
ficiently high). Nevertheless, it is useful for completeness to 
illustrate the entire range of theoretical behavior. 



Figures|28[a)-|28[c) show the fermionic density of states for 
the BEC, PG and BCS limits. These figures are important in 
establishing a precise visual picture of a "pseudogap". The 
temperatures shown are just above Tc, and for T = 0.75Tc 
and T — O.bTc. The methodology for arriving at these plots 
will be discussed in the following section, in the context of 
high Tc superconductors. Only in the BCS case is there a 
clear signature of Tc in the density of states, but the gap is so 
small and Tc is so low, that this is unlikely to be experimen- 
tally detectable. Since the fermionic gap is well established 
in the BEC case, very little temperature dependence is seen 
as the system goes from the normal to the superfluid states. 
Only the PG case, where Tc is maximal, indicates the pres- 
ence of superfluidity, not so much at Tc but once superfluid 
order is well established at T = O.STc, through the presence 
of sharper coherence features, much as seen in the cuprates. 
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FIG. 29 Typical behavior of the fermionic momentum distribution 
function at T = and Tc in the three regimes for a homogeneous 
system. The inset indicates the corresponding normalized gaps as a 
function of T /Tc. In all three regimes, Uk barely changes between 
T = Tc and T = 0. 



The inset of Fig. |29l plots the temperature dependence of 
A(T). It should be stressed that Tc is only apparent in A(r) 
in the BCS case. To underline this point, in the main body of 
Fig.|29]we plot the fermionic momentum distribution function 
rik, which is the summand in the number equation, Eq. Mi)\ , 
at T = and T = Tc. The fact that there is very little change 
from T ~ to T = Tc makes the important point that this 
momentum distribution function in a homogeneous system is 
not a good indicator of phase coherent pairing. For the PG 
case, this, in turn, derives from the fact that A(T) is nearly 
constant. For the BEC limit the excitation gap, which is dom- 
inated by /i, similarly, does not vary through Tc. In the BCS 
regime, A(r) is sufficiently small so as to be barely percepti- 
ble on the scale of the figure. 



1 . Unitary Limit and Universality 

There has been considerable emphasis on the behavior of 
fermionic gases ne ar the unitary limit where = ±oo. At- 
tention (Ho, 200l IViverit et a/.L |2004) has concentrated on 
this limit in the absence of Feshbach bosons. We now want 
to explore the more general two-channel problem here. At 
unitarity we require that 
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diverges. This will be consistent, provided 



(129) 



(130) 



It should be clear that the parameter Uq is of no particular rel- 
evance in the unitary limit. This is to be contrasted with the 
one channel problem, where Feshbach bosons are absent and 
unitarity is achieved by tuning the direct inter-fermion inter- 
action, U. 
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FIG. 30 (Color) Characteristics of the unitary limit at T = as a 
function of \go\ in the two-channel problem. Here Ubo is the den- 
sity of condensed Feshbach bosons. E'p — kp /2m defines an 
effective Fermi energy of the fermionic component with a density 
n = ri*°' — 2nbo- After rescaling, /J./ E'p, A/ E'p, and k'p^ are es- 
sentially a constant in go- This indicates that at unitarity, the system 
can be treated effectively as a two-fluid model at T = 0. The units 
for go are Ep /k^J^- 
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FIG. 31 (Color) Characteristics of the unitary limit at T = Tc as 
a function of |(7o| in the two channel problem. In contrast to the 
previous figure for T = 0, there is no rescaling which would lead to 
universal behavior of /i, A or ^ with respect to go. More than 20% of 
the system is composed of fermionic quasiparticles when \go\ —> 0. 
Here nt and Up denote the density of noncondensed Feshbach bosons 
and of the sum of FB and fermion pairs, respectively. The units for 
go are Ep/k^^. 



In the unitary Umit, the T — gap and number equations 
become 



1 + 



(131) 



and 



9o k V 



where we have used the fact that at unitarity 



^, (133) 



" [(2m-^^)C/ + .92]2 g2 

and we define Aq = Asc{T — 0). Here Uc is the critical 
coupling, given in Eq. (I17> . From these equations it follows 
that in this limit the values of /i and Ao in the ground state are 
entirely determined by the Feshbach coupling go- 

Figure |30l presents a plot of zero temperature properties 
showing their dependence on go- The two endpoints of the 
curves are of interest. For arbitrarily small but finite go the 
unitary limit corresponds roughly to /i w 0. Moreover, here 
the condensate consists exclusively of Feshbach bosons. At 
large go, (which appear to be appropriate to cuiTent exper- 
iments on ^Li and ^"K), the be havior at u nitarity is closer to 
that of the one channel problem (lHoll2004l) . Here one finds the 
so-called universal behavior with parameters ji/Ep = 0.59, 
and /S./Ef ~ 0.69 associated with the simple crossover 
ground state and a contact potential interaction^. 



To what extent is the unitary limit of the two channel prob- 
lem universal? It clearly is not in an absolute sense. However, 
there is a kind of universality implicit in these T = calcula- 
tions, as can be seen from Fig.|30]if we take into account the 
fact that there are a reduced number of fermions, as a result 
of the Bose condensate. We, thus, divide the two zero temper- 
ature energy scales by an effective Fermi energy defined by 
E'p = EF{nln^°^fl'^ with n = n*°* - 2?ig. FigurelSOlshows 
that as go decreases, the effective or rescaled Fermi energy 
E'p decreases while the Bose condensate fraction increases. 
The rescaled quantities n/E'p and A/E"^ are universal, that 
is, they are essentially independent of go, as seen in the Figure. 
A similar scaling can be done for the Cooper pair size i^. With 
a reduced effective Fermi momentum k'p = kp{n/n*°*y/'^, 
the dimensionless product k'p£^ is also independent of go. We 
can arrive at some understanding of this modified universality, 
by referring back to Sec. lV.AI In summary, the breakdown of 
universality seen in Fig.|^is to be associated with a reduced 
number of fermions, which, in turn derives from the Bose con- 
densate Once is negligible, the distinction between one 
and two channel models disappears. 

Finite temperature properties, such as Tc and A{Tc) can 
also be calculated in the unitary limit, as a function of 50- 
These are plotted in Fig.|^ Just as for the zero temperature 
energy scales, these parameters [along with /i(Tc)] increase 
monotonically with go. However, at finite temperatures there 
are noncondensed Feshbach bosons as well as noncondensed 
fermion pairs, which interact strongly with each other. These 



In our plots we find small numerical coiTections associated with our non- 



infinite cutoff. 
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FIG. 32 (Color) Characteristics of the small (;o limit as a function 
of va in the two-channel model. Unless indicated otherwise the plots 
correspond to T = 0. The molecular condensate fraction is very 
high and there is no maximum at unitarity. The gap A is very small 
(oc (?o). Here we take go ~ —lEp/k'^J^ and Uo ~ SEp/k^. 
Note the very different scales for the abscissa when comparing with 
Figs.l23landl24l 



noncondensed states are represented in the lower right inset. 
Plotted there is the fraction of noncondensed bosons [ni,{Tc)] 
and sum of both noncondensed fermion pairs and noncon- 
densed bosons [np{Tc)] at Tc- These quantities correspond 
to the expressions shown in Eqs. ( I80> and ( I92> . respectively. 

These noncondensed particles are responsible for the fact 
that at Tc a simple rescaling will not yield the universality 
found at T = 0. Moreover, these noncondensed bosons and 
fermion pairs influence calculations of T^, itself, via pseudo- 
gap effects which enter as A{Tc). These pseudogap effects, 
which are generally ignored in the literature, compete with 
superconductivity, so that Tc is bounded with its a maximum 
given by w 0.25Tf, just as seen in Fig.|53] for the one channel 
problem. Despite the fact that there is no universality, even in 
the rescaled sense, at Tc one can associate the large go limits 
with the behavior of the one channel system which has been 
rather extensively studied in the unitary limit (Ho, 2004). This 
large go limit seems to be reasonably appropriate for the cur- 
rently studied Feshbach resonances in ^Li and ''^"K. In this 
sense, these particular atomic resonances are well treated by a 
one channel problem, near the unitary limit. 

In summary, comparison between Fig.|30land FigurelsTIre- 
veals that (i) at T = 0, all Feshbach bosons and pairs are 
condensed and the two (fermionic and bosonic) condensates 
are relatively independent of each other, (ii) By contrast, 
at T = Tc, fermionic quasiparticles, noncondensed fermion 
pairs and Feshbach bosons all coexist and interact strongly 
with each other This destroys the simple scaling found at 
T = 0. 



2. Weak Feshbach Resonance Regime: I50I = 0+ 

The effects of very small yet finite go lead to a nearly 
bi-modal description of crossover phenomena. Figure |32] 
presents the behavior for very naiTow Feshbach resonances. 
It shows that there is a rather abrupt change as a function of 
magnetic field (via i^o) from a regime in which the conden- 
sate is composed only of Feshbach bosons, to one in which 
it contains only fermion pairs. Depending on the size of 
go, this transition point (taken for concreteness to be where 
2ri°/n*°* w 0.5, for example) varies somewhat, but it seems 
to occur close to the unitary limit and, importantly, for nega- 
tive scattering lengths i.e., on the BCS side. This should be 
contrasted with the large go case). This small go regime of 
parameter space was explored bv Ohashi and Griffin (Ohashi 
and Griffin, 2002). Indeed the behavior shown in the figure is 
similar to the physical picture whicl^alc^an^^too^(Falco 
and Stoof . 120041) presented. 

This figure represents a natural evolution in the behavior 
akeady indicated by Figs.l23landl24l Here, however, in con- 
trast to these previous figures the behavior of Tc is essentially 
monotonic, and the weak maximum and minimum found ear- 
lier has disappeared. For this regime, the small pseudogap 
which can be seen in the upper right panel arises predomi- 
nantly from the contribution of Feshbach bosons via Eqs. 
and ( 19 U . At the same time it should be noted that the transi- 
tion from mostly molecular bosons to mostly fermion pairs in 
the condensate is continuous. While it may look sharp par- 
ticularly for a plot with a more extended range of i^o than 
shown here, this closeup view gives a more precise picture 
of a smooth crossover transition. 

We end by noting that in this small go regime, the system is 
far from the one channel limit obtained when go is identically 
0, as shown in Fig. |25] It should be clear from the above 
discussion that Feshbach bosons play a more important role 
in the regime of unitary scattering the smaller is go, provided 
that it does not strictly vanish. 



B. Superfluidity in Traps 

Most experiments on fermionic superfluids are done in the 
presence of a trap potential V{r). For the most part, this is 
most readily included at the level of the local density approx- 
imation (LDA). In this section, we summarize the self con- 
sistent equations for spherical traps V{r) — T^uicj T . Tc is 
defined as the highest temperature at which the self consistent 
equations are satisfied precisely at the center At a tempera- 
ture T lower than Tc the superfluid region extends to a finite 
radius Rgc- The particles outside this radius are in a normal 
state. Also important is the radius i?A beyond which the ex- 
citation gap is effectively zero, to a chosen level of numerical 
accuracy. 

The self consistent equations (with simplified units h — \, 
fermion mass m = i, and Fermi energy Ep = hLu{3N)^ — 
1, where uj is the trap frequency) are given in terms of the 
Feshbach coupling constant g and inter-fermion attractive in- 
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FIG. 33 (Color) Representative spatial profiles of fermions in a harmonic trap potential at various temperatures from to T = in the three 
indicated regimes. As the system evolves from the BCS to the deep BEC regime, the fermion cloud shrinks, and a bimodal distribution of the 
density develops below Tc. Note the different scales for the vertical axes. 



teraction [/ by a gap equation 



U 



2^ - 2V{r) - V 



E 



1 - 2/(i;k) 



= . (134) 



which characterizes the excitation gap A(r) at positions r 
such that [ipair [f) — 0. Moreover, the pseudogap contribu- 
tion to A2(r) = KlJ(r) + A2^(T) is given by 

^l, = \Y.^{^'i-^Va^r) ■ (135) 



Here h represents the usual Bose function. The density of par- 
ticles at radius r: 



2 

n(j) = 2nS] + — b{ng - fiboso7i) 



+ 2J2H{1- f{Ei,))+<f{Ei,)] .(136) 



The important quantity fipair = fJ-boson is identically zero 
in the superfluid region r < Rsc, and must be solved for self 
consistently at larger radii. Here M is the bare boson mass 
equal to 2m, and the various residues, Z and Zi, and pair dis- 
persion rtq are described in Section llLCl 

The transition temperature Tc is numerically determined as 
follows: (i) A chemical potential is assigned to the center of 
the trap /i(0). The chemical potential at radius r will then be 
/i(r) — /i(0) — V{r). (ii) The gap equation [Eq. ( I134H and 
pseudogap equation [Eq. (I135H are solved at the center (set- 
ting Apg = A) to find and A(0, T^) (iii) Next the radius 
Ra is determined. (It can be defined as the point where the 
gap has a sufficiently small value, such as A(0, Tc)/1000). 
(iv) Following this, the pseudogap equation (Eq. (I135H is 
solved for A(r, Tc) up to the radius Ra- Then n{r) is deter- 
mined using Eq. (I136> . (v) Finally, n{r) is integrated over all 
space ( J cPrn{r)) and compared with the given total number 
of fermions N. The procedure is repeated until the prescribed 
accuracy has been reached. To extend this algorithm below Tc 
a similar approach is used which now includes an additional 
step in which Rsc is computed. 



Figure|33]shows a plot of the density profiles obtained from 
solving the self consistent equations. Here the Thomas-Fermi 
radius, called Rtf, is calculated for the non-interacting case. 
In the BCS regime there is virtually no temperature depen- 
dence in the density profile as expected from the small val- 
ues for A. In the unitary and BEC regimes as the system ap- 
proaches Tc, n{r) can be reasonably well fit to a Gaussian. A 
bi-modality in the distribution is evident in the BEC regime, 
thus providing a clear signature of the superfluid phase. How- 
ever, in the PG or unitary it is more difficult to establish super- 
fluidity, since the profiles do not show a clear break at the po- 
sition where the condensate first appears (Stajic et ai, 2005ai). 

These observations differ some what from other theoreti- 
cal pr edictions in the literature dChiofalo et al\. l2002l iHol 
l2004h . The differences between the present theory and this 
other work can be associated with the inclusion of noncon- 
densed pairs or pseudogap effects. These noncondensed pairs 
smooth out the otherwise abrupt transition between the con- 
densat e and the fermionic excitations. In contrast to earlier 
work (iPerali et aL. .,2004b,) these st udies, wh ich are based on 
the standard crossover ground state ( lLeggetd. fT980). yield only 
monotonic density distributions. Future experiments will be 
needed to look f or the secondary max ima predicted by Strinati 
and co-workers dPerali et al.l l2004ah in the BEC and unitary 
regimes. 

Finally, Fig.|34]shows the behavior of Tc (inset) in a trap, 
together with the superfluid density (main figure). The latter 
is computed following the discussion in Sec. [Ill] From the 
inset one can rea d off that T^jEp » 0.3 at resonance, which 
is very close to (Perali et ai, 2004a), or slightly smaller than 
other estimates (.Ohashi and Griffin, .2003) in the literature. 

The structure seen in the plot of Tc for the trapped gas re- 
flects the minimum and maximum found for the homogeneous 
solution (Fig. I23> : these are related, respectively, to where 
/i w and where unitary scattering occurs. The Tc curve fol- 
lows generally the decrease in the density at the center of the 
trap from BEC to BCS. The temperature dependences of the 
(integrated over r) superfluid density reflect the nature of ex- 
citations of the condensate. The BEC limit has only bosonic, 
while the BCS has only fermionic excitations. The unitary 
case is somewhere in between. The fermionic contributions 
are apparent in the BCS regime, where Ng decreases almost 
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FIG. 34 Behavior of superfluid density A'^s as a function of T in the 
three different regimes. Note that the negative curvatures at low T 
for the unitary and BEC cases are different from the linear behavior 
in the BCS case. In the BEC limit, - iV^ oc T^/^. In the inset 
we show Tc as a function of magnetic detuning vq in a tras. We 
have taken Uq = —O.dEp/k'p and go — —SSEp/k^^, ba^ed on 
experimental data for ''"K. 



linearly at low T. For an s-wave superfluid this refledts the 
contributions from low A(r) < T regions, near the edge of 
the trap. These are present to a lesser degree in the unitary 
case. The bosonic contribution to Ng and to thermodyii amies 
(in both the BEC and unitary regimes) is governed by the im- 
portant constraint that Upairii") vanishes at positions r where 
there is a co-existent superfluid. A similar constraint is im- 
posed at the level of Har tree-Fock-Bogoliubov theory for true 
bosons dShi and Grif fin. 1998). These gapless pair excitations 
are responsible for a T^/'^ temperature dependence in Ng. In 
this way the low T power law dependences are similar to their 
analogues in the homogeneous case of Section rVl.DI m ich as 
found for true Rose gases ( Carr et al. .2004) in a trap, i 

As for the BEC regime there is an interesting obserl/ation 
to make about the contribution to the density profile associ- 
ated with the condensate. Experimentalists generally assume 
a Gaussian form throughout the trap for the nonconcensed 
pairs. In this wa y they estimate the condensate fraction Here 
we find dStaiic et al., 2005a) that this Gaussian form is not 
correct in the superfluid region of the trap, and that this Gaus- 
sian approximation will underestimate (by w a factor c f 1.5) 
the condensate fraction in the BEC regime. 



C. Comparison Between Theory and Experiment: 
Atoms 



As time progresses it will be i mportant to as sess 
sumed ground state wavefunction (Leggett, '1980') and 
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channel va riant (Milstein et al., 2002; Ohashi and Griffin, 
|2002; Stai ic~^a/.l l2005dh in more detail. A first gijnera- 
tion set of measurements and theories have indicated t n ini- 
tial level of success in addressing four experiments: RI' pair- 
ing gap spectroscopy and collective mode measurements, as 
well as more recently, the particle density profiles dStaiic 



and thermodynamics ( iKinast et al\. l2005h of 
Fermf gas es in the unitary r egime. The first of these ex- 
ents dChin et oil 120041) seem to provide some support 
for tie fact that the Bogoliubov excitations with dispersion, 
i^k = \/ (ck — + A2(T) have a "gap parameter" or pseu- 
dogaj), A(T) which is in general different from the supercon- 
ducting order parameter, as a consequence of noncondensed 
pairs at and above Tc- Torma and colleagues ( Kin nunen et aZj 
1200 41:)) have used the formalism outlined i n Secti on III. Bl to 
address RF spectroscopy. Subsequently (^ Kinnunen et al\ 
5004;i) they have combined this analysis with an approximate 
LDA t reatment of the tra p to make contact with RF experi- 
ment! (IChin et al l l2004h . Important to their work is the fact 
that tlie edges of a trap have a free atom character, where there 
is a vanishing A(r). 

Additional support for the wavefunction comes from the 
collective mode behavior near T — Q. Th e mode frequen - 
cies can be connected to an equation of state ( Stringari, 2004). 
Numerical studies ( Heiselberg. 2004b; Hu et aL 2004) of this 
equal ion for the standard crossover ground state demonstrate a 
good fit to the measured breathing mode frequencies (Barten- 
stein et al., 2004a; Kinast et al., 2004). Importantly, in the 
near- iEC limit the frequencies decrease with increasing mag- 
netic field; this is opposite to the behavior predicted for true 
bosoi s (Stringari, 2004), and appears to reinforce the notion 
that for the magnetic fields studied thus far, the "BEC" regime 
should be associated with composite rather than real bosons. 

At iT = 0, calculations of the collective mode frequencies 
are based on an equation of state. Moreover, the zero temper- 
ature equation of state associated with the mean field ground 
state :an be computed analytically in the near-BEC limit 
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where the first two terms on the right hand side were pre- 
viously discussed following Eq. (I55> . Importantly, the third 
term is responsible for the observed agreement between the- 
ory and collective mode experiments, in the near-BEC regime. 

A serious shortcoming of the wavefunction ansatz has been 
emphasized in the literature. This is associated with the ra- 
tio ol the inter-boson to the fermionic scattering length near 
the unitary regime, but on the BEC side of the Feshbach reso- 
nanctt. Experimen t ally the numbers are considera bly smaller 
llBart^nstein et a/.L l2004hl: Iz ^ierlein et fl/.L l2004h (0.6) than 
the fector of 2 predicted for the one channel model discussed 
belov Eq. \55l . They a re consistent, howeve r, with more exact 
few body calculations llPetrov et a/.L l2004l) . These few body 
calculations apply to a regime of magnetic fields where kpcis 
is smaller than about 0.3 -0.4. This seems to be close to the 
bouniiary for th e "quasi-BEC" regime wh ich can be accessed 
expeiimentally ( Bartenstein et a/.Ll2004bh in ^Li. 

It is interesting to repeat the T = many body calcula- 
tions now with the inclusion of FB, since these two channel 
effects begin to be important at kpas smaller than around 0.2 
-0.25. Here, one sees (Staiic et al, 2005b), that the factor of 2 
found in the one channel calculations is no longer applicable, 
and th e values of the ratio (of bosonic to fermionic scattering 
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FIG. 35 Phase diagram for a quasi-two-dimensional (2D) d-wave 
superconductor on a lattice. Here the horizontal axis corresponds to 
the strength of attractive interaction —U/it, where t is the in-plane 
hopping matrix element. 



lengths) are considerably smaller, and, thus, tending in the di- 
rection which is more in line with experiment. This reflects 
the decrease in the number of fermionic states which mediate 
the interaction between molecular bosons. 

In this way the mean field theory presented here looks 
promising. Nevertheless, uncertaint ies remain. Current theo- 
retical work jKinnunen et a/.LE004ah on the RF experiments is 
limited to the intermediate coupling phase, while experiments 
are also available for the near-BEC and near-BCS regimes, 
and should be addressed theoretically. In add ition, the pre- 
dictions of alternative approaches ( Pera li et al. l|2004b) which 
posit a different ground state will need to be compared with 
these RF data. The col lective mode calculations ( Heiselbera, 
l2004hl:lHu et a/.L Eo04l) which have been performed have not 
yet considered the effects of Feshbach bosons, which may 
change the behavior in the near-BEC regime. 

It is, finally, quite possible that incomplete T — Q conden- 
sation will become evident in future experiments. If so, an al- 
ternative wavefunc tion will have to be c ontemplated. Indeed, 
it has been argued llHollandefaa i2004) that the failure of the 
one chan nel model to match the calculations of Petrov and co- 
workers (iPetrov et aU |2p04) is evidence that there are serious 
weaknesses. At an experimental level, new pairing gap spec- 
troscopies aoDear^^^merginj^^^^M^ragi^ac^Bruun 
and Baym, 2004 lGreiner et a/.Ll2005h which will further test 
these and subsequent theories. 



VI. PHYSICAL IMPLICATIONS: HIGH 
SUPERCONDUCTIVITY 

We begin with an important caveat that many of the exper- 
iments outlined in this section can a lternatively be addressed 
from a Mott physics point of view ( iLee et ad 12004). Here 
we summarize the interpretation of these experiments within 
BCS-BEC crossover theory with the hope that this discussion 
will be of value to the cold atom community, and that ulti- 
mately, it will help to establish precisely how relevant is the 
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FIG. 36 Inverse pair mass (solid line) and Tc (dashed line) as a func- 
tion of coupling strength — (7/4t for a d-wave superconductor on a 
quasi-2D square lattice, where t is the in-plane hopping matrix ele- 
ment. Tc vanishes when the pair mass diverges, before the bosonic 
regime is reached. The electron density is taken to be n = 0.85 per 
unit cell, corresponding to optimal hole doping in the cuprate super- 
conductors. 



BCS-BEC crossover scenario for understanding high temper- 
ature superconductivity. 



A. Phase Diagram and Superconductor-Insulator 
Transition: Boson Localization Effects 

The high Tc superconductors are different from the ultra- 
cold fermionic superfluids in one key respect; they are d-wave 
superconductors and their electronic dispersion is associated 
with a quasi-two dimensional tight binding lattice. In many 
ways this is not a profound difference from the perspective of 
BCS-BEC crossover. Figure l35l shows a plot of the two im- 
portant temperatures Tc and T* as a function of increasing at- 
tractive coupling. On the left is BCS and the right is PG. The 
BEC regime is not apparent. This is because Tc disappears 
before it can be accessed. At the point where Tc vanishes, 
^/Ef « 0.8. 

A competition between increasing T* and Tc is also appar- 
ent in Fig. 1^ This is a consequence of pseudogap effects. 
There are fewer low energy fermions around to pair, as T* 
increases. It is interesting to compare Fig. |35l with the exper- 
imental phase diagram plotted as a function of x in Fig.0 If 
one inverts the horizontal axis (and ignores the AFM region, 
which is not addressed here) the two are very similar To make 
an association between the coupling U and the variable x, it 
is reasonable to simply fit T*{x). We do not dwell on this last 
step here, since we wish to emphasize crossover effects which 
are not complicated by "Mott physics". 

Because of quasi-two dimensionality, the energy scales of 
the vertical axis in Fig.|35]are considerably smaller than their 
three dimensional analogues. Thus, pseudogap effects are 
intensified, just as conventional fluctuation effects are more 
apparent in low dimensional systems. This may be one of 
the reasons why the cuprates are among the first materials to 
clearly reveal pseudogap physics. Moreover, the present cal- 
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FIG. 37 Evolution of spectral functions of a pseudogapped superconductor with temperature. As soon as T drops below Tc, the spectral peaks 
split up rapidly, signalling the onset of phase coherence. 



culations show that in a strictly 2D material, Tc is driven to 
zero, by bosonic or fluctuation effects. This is a direct reflec- 
tion of the fact that there is no Bose condensation in 2D. 

Figure |36] presents a plot (Chen et al., 1999) of boson lo- 
calization effects associated with d-wave pairing. A diver- 
gence in the effective mass m/AI* occurs when the cou- 
pling strength, or equivalently T* becomes sufficiently large. 
Once the bosons localize the system can no longer support 
superconductivity and Tc vanishes. It is generally believed 
that a system of bosons with short range repulsion is either 
a superfluid or a Mott insulator in its ground state (Fisher 
et fl/., 1989). Moreover, the insulating phase is presumed to 
contain some type of long range orde r. Interes t ingly, there 
are recent reports of ordered phases (Fu et al.\ '20041 Ver- 
shinin et al., 2004) in the very underdoped cuprates. These 
appear on the insulating side and possibly penetrate across 
the superconductor-insulator boundary; they hav e been inter- 
preted by some as Cooper pair dens i ty waves (|Chen et al.l 
|2004; Me Ukvan and Tesanovi c. 2004; Tesanovic, 2004|). The 
observation of boson localization appears to be a natural con- 
sequence of BCS-BEC theory; it should be stressed that it is 
associated with the d-wave lattice case. Because d-wave pairs 
are more spatially extended (than their s- wave counterparts 
they experience stronger Pauli repulsion (IChen et al.l Il99 
and this inhibits their hopping. 



B. Superconducting Coherence Effects 

The presence of pseudogap effects raises an interesting set 
of issues suiTounding the signatures of the transition which 
the high Tc community has wrestled with, much as the cold 
atom community is doing today. For a charged superconduc- 
tor there is no difficulty in measuring the superfluid density, 
through the electrodynamic response. Thus one knows with 
certainty where Tc is. Nevertheless, people have been con- 
cerned about precisely how the onset of phase coherence is 
reflected in thermodynamics, such as Cv or in the fermionic 
spectral function. One understands how phase coherence 
shows up in BCS theory, since the ordered state is always ac- 
companied by the appearance of an excitation gap. When a 
gap is already well developed at Tc, how do these signatures 



emerge 

To address these coherence effe cts one has to intr oduce a 
distinction between the self energy iChen et a/.L l2001') associ- 
ated with noncondensed and condensed pairs. This distinction 
is blurred by the approximation of Eq. ( I89> . Above, but near 
Tc, or at any temperature below we now use an improved ap- 
proximation (Malv et al. . 1999a b) 
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This is to be distinguished from Yigc where the condensed 
pairs are infinitely long-lived and there is no counterpart for 
7. The value of this parameter, and even its T-dependence is 
not particularly important, as long as it is nonzero. 

Figurel^plots the fermionic spectral function at = 0, 
called A{uj), as the system passes from above to below Tc. 
One can see in this figure that just below Tc, A{uj) is zero 
at a point uj = 0, and that as temperature further decreases 
the spectral function evolves smoothly into approximately two 
slightly broadened delta functions, which are just like their 
counterparts in BCS. In this way there is a clear signature as- 
sociated with superconducting coherence. To compare with 
experimental data is somewha t complicated, since me asure- 
ments of the spectral function ( Campuzano et al.\*2004\ Dam- 
ascelli et al., 2003) in the cuprates also reveal other higher 
energy features ("dip, hump and kink"), not specifically asso- 
ciated with the effects of phase coherence. Nevertheless, this 
Figure, like its experimental counterpart, illustrates that sharp 
gap features can be seen in the spectral function, but only be- 
low Tc. 

Analogous plots of superconducting coherence effects are 
presented in Fig. |38]in the context of more direct compari- 
son with experiment. Shown are the results of specific heat 
and tunneUng calculations and their experimental counterparts 
("Renn er et al 1 119981: iTallon and Lorami 1200 ll) . The latter 
measures, effectively, the density of fermionic states. Here 
the label "PG" coiTesponds to an extrapolated normal state 
in which we set the superconducting order parameter Age to 
zero, but maintain the the total excitation gap A to be the same 
as in a phase coherent, superconducting state. Agreement be- 
tween theory and experiment is satisfactory. We present this 
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FIG. 38 Superconducting state (SC) and extrapolated normal state (PG) contributions to SIN tunneling characteristics dl /dV (upper panels) 
and thermody namics S/T and Cy /T (lower panels). We compare theoretical calculations (left) with experiments (right) on tunneling for 
BSCCO from <Renner ef 0/1119981) and on specific heat for Yo.8Cao.2Ba2Cu307-« (YBCO) from jTallon and Loranll200ll) . The theoretical 
SIN curve is calculated for T = Tc/2, while the experimental curves are measured outside (dashed line) and inside (solid line) a vortex core. 
The dashed and solid curves in the upper left panel are for SC and PG state, respectively. In both lower panels, the solid and dashed curves 
indicate the SC and PG state, respectively. 
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FIG. 39 Comparison of normalized superfluid density ps / ps (0) ver- 
sus T/Tc for experiment (main figure) and theory (inset) in high Tc 
superconductors with variable doping concentration x. Both exper- 
iment and theory show a quasi-universal behavior with respect to 
doping. 
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FIG. 40 Offset plot of superfluid density ps comparing experiment 
(main figure) with theory (inset). Both experiment and theory show 
the same quasi-universal slope dps/dT for different hole concentra- 
tions X. 



artificial normal state extrapolation in discussing thermody- 
namics in order to make contact with its experimental coun- 
terpart. However, it should be stressed that in zero magnetic 
field, there is no coexistent non-superconducting phase. BCS 
theory is a rather special case in which there are two possible 
phases below and one can, thereby, use this coexistence to 
make a reasonable estimate of the condensation energy. When 
Bose condensation needs to be accommodated, there seems to 
be no alternative "normal" phase below T,- 



C. Electrodynamics and Thermal Conductivity 

In some ways the subtleties of phase coherent pairing are 
even more perplexing in the context of electrodynamics. Fig- 
ure El presents a paradox in which the excitation gap for 
fermions appears to have little to do wit h the behavior of the 
superfluid density. To address these data dStaiic et al l l2003bh 
one may use the formalism of Sec. 1111. Al One has to introduce 
the variable x (which accounts for Mott physics) and this is 
done via a fit to T*{x) in the phase diagram. Here for Fig. 140! 
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it is also necessary to fit ps{T = 0, x) to experiment, although 
this is not important in Fig. 03 The figures show a reasonable 
correspondence (Stajic et ai, 2OO2B) with experiment. The 
paradox raised by Fig.^jis resolved by noting that there are 
bosonic excitations of the condensate, as in Eq. ( II 13> and that 
these become more marked with underdoping, as pseudogap 
effects increase. In this way ps does not exclusively reflect the 
fermionic gap, but rather vanishes "prematurely" before this 
gap is zero, as a result of pair excitations. 

The optical conductivitv a(uj) is similarlv modified (Iyen- 
gar et ai, 2003; Timusk and Statt, 1999). Indeed, there is an 
intimate relation between ps and (t{u!) known as the f-sum 
rule. Optical conductivity studies in the literature, both theory 
and experiment, have concentrated on the low uj,T regime 
and the interplay between impurity scattering and d wave su- 
perconductivity ( Berlin s kv et ai. 2000). Als o of interest are 
unus ually high ui tails jMoleeraaf et adl2002l Santander-Svro 
et flZ.7'2003') in the real part of (j{uj) which can be inferred 
from sum-rule arguments and experiment. Figure Elraises a 
third set of questions which pertain to the more global behav- 
ior of cr. In the strong PG regime, where A has virtually no T 
dependence below Tc, the BCS-computed a{uj) will be sim- 
ilarly T-independent. This is in contrast to what is observed 
experimentally where a{uj) reflects the same T dependence as 
in ps{T), as dictated by the f-sum rule. 

One may deduce a consequence of this sum rule, based on 
Eq. ( II 13> . If we associate the fermionic contributions with the 
first term in square brackets in this equation, and the bosonic 
contributions with the second. Note that the fermionic trans- 
port terms reflect the full A just as do the single particle prop- 
erties. (This may be seen by combining the superconducting 
and Maki Thompson or PG diagrams in Appendix IbI. We 
may infer a sum rule constraint on the bosonic contributions, 
which vanish as expected in the BCS regime where Apg van- 
ishes. We write 

2 f°° , A^ /T7 \Bcs 

- df7(T^°^°"''(f7,T) = ( — ) (T). (138) 

TT 7g A^ V TO / 

The bosonic contributions can be determined most readily 
from a framework such as time dependent Ginzburg Landau 
theory, which represents rather well the contributions from the 
Aslamazov-Larkin diagrams, discussed in Appendix 151 The 
bosons make a maximum contributio n at Tr,. 

The resulting optical conductivity ("ivengar et all l2003h in 
the superconducting state is plotted in Fig. l41l below. The key 
feature to note is the decrease in width of the so-called Drude 
peak with decreasing temperature. In the present picture this 
pseudogap effect can be understood as coming from the con- 
tribution of bosonic excitations of the condensate, which are 
present at Tc, but which disappear at low T. This behavior 
is reasonably consistent with what is reported experimentally 
dTimuskand Stattlll999l) . where "the effect of the opening of 
a pseudogap is a narrowing of the coherent Drude peak at low 
frequency". The calculations leading to Fig.|^ indicate that 
by 20 K most of the conductivity is fermionic except for a 
high frequency tail associated with the bosons. This tail may 
be responsible in part for the anomalously hig h uj contribu- 
tions to a required to satisfy the sum rule ( Molegraaf et all 
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FIG. 41 Real part of ac conductivity at low and high T below Tc. 
Here pseudogap effects, which insure that A(r) ~ constant below 
Tc, require that on e go beyond BC S theory to explain the data shown 
in the inset 1 Puchkov et all\T99(il . In the theory, inelastic and other 
effects which presumably lead to the large constant background term 
observed in the data have been ignored. 

'2002^.'Sant ander-Syro et adl2003[) . 

Although the normal state calculations have not been done 
in detail, once the temperature is less than T*, pseudogap ef- 
fects associated with the fermionic contribution will also lead 
to a progressive narrowing in the Drude peak with decreasing 
T, down to Tc, due to the opening of an e xcitation gap. This 
mechanism (which is widely ascribed to dTimusk and Stattj 
1 1 9991) in the literature) reflects longer lived electronic states, 
associated with the opening of this gap. However, it cannot be 
applied to explain the progressive narrowing of the conductiv- 
ity peaks, once A(T) becomes temperature independent, as is 
the case at low T. In this way, one can see that pseudogap ef- 
fects show up in the optical conductivity of underdoped sam- 
ples, both from the fermionic contributions (near and below 
T*) and the bosonic contributions (below Tc) as a narrowing 
of the Drude peak. 

The data shown in the inset of Fig.|^reveal that the con- 
densate is made from frequency contributions primarily up to 
w 0.15 eV. These data pose a challenge for applying strict 
BCS theory to this problem, since one consequence of a tem- 
perature independent A in the underdoped phase is the ab- 
sence of T dependence in the optical conductivity. This is the 
same challenge as posed by Fig.^] Strikingly there is no ex- 
plicit "gap feature" in the optical conductivity in the normal 
state. Such a gap feature, that is a near vanishing of a for 
an extended range of frequencies, will occur in models where 
the pseudogap is unrelated to superconductivity; it is not seen 
experimentally, thus far 

By contrast with the electrical conductivity, the thermal 
conductivity is dominated by the fermionic contributions at 
essentially all T. This is because the bosonic contribution 
to the heat current (as in standard TDGL theory (Larkin and 
Varlamov, 2001)), is negligibly small, reflecting the low en- 
ergy scales of the boson s. Thermal conductivity, experiments 
( [Sutherland et all l2003h in the high Tc superconductors pro- 
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vide some of the best evidence for the presence of fermionic 
d-wave quasiparticles below T^,. In contrast t o the ac con duc- 
tivity, here one sees a universal low T limit 1 19931) . and 
there is little to suggest that something other than conventional 
d-wave BCS physics is going on here. This cannot quite be the 
case however, since in the pseudogap regime, the temperature 
dependence of the fermionic excitation gap is highly anoma- 
lous, as shown in Fig.|3] compared to the BCS analogue. 



D. Thermodynamics and Pair-breaking Effects 

The existence of noncondensed pair states below Tc affects 
thermodynamics, in the same way that electrodynamics is af- 
fected^_as_discu^se^above^,^orw 

et a/.. l2000l) that the dispersion will lead to ideal Rose 
gas power laws in thermodynamical and transport properties. 
These will be present in addition to the usual power laws or 
(for s-wave) exponential temperature dependences associated 
with the fermionic quasiparticles. Note that the dependence 
reflects the spatial extent ^pg, of the composite pairs , and this 
size effect has no natural counterpart in true Rose systems. It 
should be stressed that numerical calculations show that these 
pair masses, as well as the residue Z are roughly T indepen- 
dent constants at low T. As a result, Eq. \9A\ implies that 

A2^ = A2(r)-AL(T)(xr3/2. 

The consequence s of these observations can now be listed 
dChen et al\. l2000l) . For a quasi-two dimensional system, 
Cy/T will appear roughly constant at the lowest temperatures, 
although it vanishes strictly at T = as T^/^. The superfluid 
density ps (T) will acquire a T'^/^ contribution in addition to 
the usual fermionic terms. Ry contrast, for spin singlet states, 
there is no explicit pair contribution to the Knight shift. In 
this way the low T Knight shift reflects only the fermions and 
exhibits a scaling with T/A(0) at low temperatures. Experi- 
mentally, in the cuprates, one usually sees a finite low T con- 
tribution to Cv/T. A Knight shift scaling is seen. Finally, also 
observed is a deviation from the predicted d-wave linear in T 
power law in ps- The new power laws in Cy and ps are con- 
ventionally attributed to impurity effects, where ps cx T^, and 
Cv/T (X const. Experiments are not yet at a stage to clearly 
distinguish between these two alternative explanations. 

Pair-breaking effects are viewed as providing important in- 
sight into the origin of the cuprate pseudogap. Indeed, the 
different pair-breaking sensitivities of T* and Tc are usually 
proposed to support the notion that the pseudogap has nothing 
to do with superconductivity. To counter this inference, a de- 
tailed set of studies was conducted, (based on the REC-RCS 
scenario), of pair-br eaking in the presence of impurities ( Chen 
and Schrie ffer, 2002 HKao et aZ.Ll2002h and of magnetic fields 
( iKao et all i2001i) . These studies make it clear that the su- 
perconducting coherence temperature Tc is far more sensitive 
to pair-breaking than is the pseudogap onset temperature T* . 
Indeed, the phase diagram of Fig. |35l which mirrors its experi- 
mental counterpart, shows the very different, even competing 
nature of T* and Tc, despite the fact that both arise from the 
same pairing correlations. 



E. Anomalous Normal State Transport: Nernst and Other 
Precursor Transport Contributions 

Much attention is given to the anomalous behavior of the 
Nernst coefficient in the cuprates (Xu et ai, 2000). This coef- 
ficient is rather simply related to the transverse thermoelectric 
coefficient a^y which is plotted in Fig. ^] In large part, the 
origin of the excitement in the literature stems from the fact 
that the Nernst coefficient behaves smoothly through the su- 
perconducting transition. Relow Tc it is understood to be as- 
sociated with superconducting vortices. Above Tc if the sys- 
tem were a Fermi liquid, there are arguments to prove that 
the Nernst coefficient should be essentially zero. Hence the 
observation of a non-negligible Nernst contribution has led to 
the picture of normal state vortices. 

Another way to view these data is that there are bosonic 
degrees of freedom present in the normal state, as in the RCS- 
REC crossover approach. To describe this dynamics we begin 
with the TDGL scheme outlined in Sec. lI.H.ll This formalism 
must be modified somewhat however When "pre-formed" 
pairs are present at temperatures T* high compared to Tc, 
the classical bosonic fields of convention al TDGL must be re - 
placed by their quantum counterparts (T an and LevinlE'004^ ■ 
We investigate this fairly straightforward extension of TDGL 
in the next sub-section, and subsequently explore the resulting 
imphcations for transport. 



1 . Quantum Extension of TDGL 

To make progress on the quantum analogue of TDGL we 
study a Hamiltonian describing bosons coupled to a parti- 
cle reservoir of fermion pairs. The latter are addressed in 
a quantum mechanical fashion. Our treatment of the reser- 
voir has strong similarities to the approach of Caldeira and 
Leggett ( Caldeira and Legaett, 1983). These bosons are in the 
presence of an electromagnetic field, which interacts with the 
fermion pairs of the reservoir as well, since they have charge 
e* = 2e. For simplicity we assume that the fermions are dis- 
persionless. This Hamiltonian is given by 

H = ^£imV'/(i)exp(-ze*am(i))V'm(t) 

Im 

1 ii 

+ J2{{b, - e*cl))vlv, + Qij^vj + O.^}. (139) 
ii 

Here £im is the hopping matrix element for the bosons. 
Cim(^) is the usual phase factor associated with the vector 
potential. Annihilation operators for the reservoir, Wi and 
Vi (with infinitesimal coupling constants rji and Q), represent 
positive and negative frequencies respectively and need to be 
treated separately. The energies and bi of the reservoir are 
both positive. 

The bosonic propagator T(k, oj) can be exactly computed 
from the equations of motion. Here A(k, oj) — Re2zT(k, uj) 
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FIG. 42 Calculated transverse thermoelectric response, which ap- 
pears in the Nemst coefficient. Inset plots experimental iXu et all 
l2000t) data. Both the theoretical and experimental values of a^y are 
normalized by the value of a^y for an optimally doped sample at 
2n, 0.0232(B/T) V/(Knm). 
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FIG. 43 dc resistivity p (normalized by pg = s/e^) vs. T, for 
quasi-2D underdoped (UD) samples with two different forms for 
p-pmr. These are shown plotted in the left inset as = — /ipair/Im7. 
The right inset shows corresponding axy, and arrows indicate Tc- 



is the boson spectral function, and 

T(k,c^)= (cc;-e(k)-I]i(cc;) + ^I]2(c^))~' (140) 
We have 

E^iu) = 27r^|77,|25(c^-a,)-27r^|GP'^(c^ + 6.), (141) 

i i 

with Si (ti)) defined by a Kramers Kronig transform. 

The reservoir parameters a^, hi, rji and Q which are of no 
particular interest, are all subsumed into the boson self energy 
E2(w). From this point forward we may ignore these quan- 
tities in favor of the boson self energy. We reiterate that this 
theory makes an important simplification, that the reservoir 
pairs have no dispersion. For this case one can solve for the 
exact transport coefficients. 

For small, but constant magnetic and electric fields and 
thermal gradients we obtain the linearized response. For no- 
tational convenience, we define tiab(k) = g^'^gfe Then 
we can write a few of the in-plane transport coefficients which 
appear in Eqs. (|8} and (|9} as 

c*2 



a 



d^fc dw ,,,, .db(uj) 



(142) 



and 



= *2 1 



67 J (27r)3 2tt 



VxVxVyyUjA^{'k,Lu) 



db{uj) 



(143) 



These equations naturally correspond to their TDGL coun- 
terparts (except for different phenomenological coefficients) 
when T w T^. 



within the framework of BCS-BEC crossover. The results for 
a^y, which relates to the important Nernst effect, are plotted 
in Fig. 1^ with a subset of the data plotted in the upper right 
inset. It can be seen that the agreement is reasonable. In this 
way a "pre-formed pair" picture is, at present, a viable alter- 
native to "normal state vortices". 

Within the crossover scenario, just as in TDGL, there are 
strong constraints on other precursor transport effects: para- 
conductivity, diamagnetism and optical conductivity are all in- 
directly or directly connected to the Nernst coefficient, in the 
sense that they all derive from the same dynamical equation of 
motion for the bosons. An alogous inferences have been made 
jOeshkenbein et al\ . \l99lh . based on similar bosonic general- 
izations of TDGL. 

The calculated b ehavior of a number of these other proper- 
ties is summarized ( iTan and "Le vin. 2004) in Figs.l43land l44l 
below, which indicate both quasi-two dimensional and more 
three dimensional behavior. In the plots for the resistivity, 
p, a fairly generic (Leridon et ai, 2001) hnear-in-7 contri- 
bution ^ from the fermions has been added, as a background 
term. A central theme of concern in the literature is whether 
one can understand the behavior of p simultaneously with that 
of a.r „. Experiments jLeridon et all l20()ll IWatanabe et all 
Il997^ indicate that there is a small feature corresponding to 
7* at which p deviates from the background linear fermionic 
term, but that the deviation is not noticeably large except in 
the immediate vicinity of 7c. This is in contrast to the behav- 
ior of axy which shows a significant high temperature onset. 
A key observation from the theoretical calculations shown in 
the figures is that the onset temperatures for the Nernst coef- 
ficient and the resistivity are rather different with the former 
being considerably higher than the latter Moreover all onsets 
appear to be significantly less than 7* . Qualitative agreement 



2. Comparison with Transport Data 

The formalism of the previous sub-section and in particu- 
lar, Eqs. il42\ and ( I143> . can be used to address transport data 



' In reality, one should include an additional temperature dependence from 
pseudogapped fermions which presumably has a small upturn, relative to a 
linear curve, as T decreases. 
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FIG. 44 dc resistivity and a^y, for an underdoped, more thiree- 
dimensional, system like YBCO. Tlie arrow indicates Tc. The mag- 
netization M is plotted as the inset, and is independent of dimension- 
ality. 

with experiment is not unreasonable. However, a more quan- 
titative comparison with experiment has not been established, 
since it depends in detail on a phenomenological input: the 
temperature dependence of the pair chemical potential /ipaii . 

Another important feature of these plots is the fact that there 
is virtually no diamagnetic precursor, as was also found to be 
the case in simple TDGL. To see this, note that the magne- 
tization can be written LSchmid. .1969: .Tan and Levin. .20041) 

B 3 he Vl + (2^c/s)2' 

where a is the fine structure constant, c is the speed of light, s 
is the interlayer spacing, and ^afc (Cc) are coherence lengths in 
the direction parallel (perpendicular) to the planes. From this 
result we can estimate the order of magnitude of the diamag- 
netic susceptibility due to preformed pairs. For a typical high 
Tc superconductor, T ^ 50K, so that A = 3 x lO^^'m, 

and A/a ~ 4 X IG^^m. Because the coherence lengths and 
the interlayer spacing are on the nanometer length scale, pre- 
cursor effects in the diamagnetic susceptibility are of the order 
10^^, and thus extremely small, as is consistent with experi- 
ment ( Wang etal. 2001). 

In order to make the BCS-BEC scenario more convincing 
it will be necessary to take these transport calculations below 
Tc- This is a project for future research and in this context it 
will ultimately be important to establish in this picture how 
superconducting state vortices are affected by the persistent 
pseudogap. Magnetic field ind uced magne t ic ord er which has 
been associated with vortices llLake et al\ |2002) will, more- 
over, need to be addressed. The analogous interplay of vor- 
tices and pseudogap will also be of interest in the neutral su- 
perfluids. 

VII. CONCLUSIONS 

In this Review we have summarized a large body of work 
on the subject of the BCS-BEC crossover scenario. In this 



context, we have explored the intersection of two very differ- 
ent fields: high superconductivity and cold atom superflu- 
idity. Theories of cuprate superconductivity can be crudely 
classified as focusing on "Mott physics" which reflects the 
anomalously small zero temperature superfluid density and 
"crossover physics", which reflects the anomalously short co- 
herence length. Both schools are currently very interested 
in explaining the origin of the mysterious pseudogap phase. 
In this Review we have presented a case for its origin in 
crossover physics. The pseudogap in the normal state can 
be associated with meta-stable pairs of fermions; a (pseudo- 
gap) energy must be supplied to break these pairs apart into 
their separate components. The pseudogap also persists be- 
low Tc in the sense that there are non condensed fermion pair 
excitations of the condensate. These concepts have a natural 
analogue in self consistent theories of superconducting fluctu- 
ations, but for the crossover problem the width of the "critical 
region" is extremely large. This reflects the much stronger- 
than-BCS attractive interaction. 

It was not our intent to shortchange the role of Mott physics 
which will obviously be of importance in our ultimate under- 
standing of the superconducting cuprates. There is, however, 
much in this regard which is still uncertain associated with 
establishing the simultaneous relevance and existence of spin- 
charge separation (Anderson, 1997 ), stripes ([K ivelson et all 
120031) . and hidden order parameters IChakravartv et a/.Ll200ll) . 
What we do have in hand, though, is a very clear experimen- 
tal picture of an extremely unusual superconductor in which 
superconductivity seems to evolve gradually from above Tc to 
below. We have in this Review tried to emphasize the common 
ground between high Tc superconductors and ultracold super- 
fluids. These Mott issues may nevertheless, set the agenda for 
future cold atom sUidies of fermions in optical lattices ( Hof- 
stetter et ai, I2OO2I) . 

The recent discovery of superfluidity in cold fermion gases 
opens the door to a new set of fascinating problems in con- 
densed matter physics. Unlike the bosonic system, there is 
no ready-made counterpart of Gross-Pitaevskii theory. A new 
mean field theory which goes beyond BCS and encompasses 
BEC in some form or another will have to be developed in 
concert with experiment. As of this writing, there are four 
experiments where the simple mean field theory discussed in 
this review is in reasonable agreement with the data. These 
are the collective mo de studies over the entire range of ac ces- 
sible magnetic fields (lHeiselbergir2004bl: lHu et alW 2004 . In 
addition in the unitary regime, RF spectroscopy-based pairing 
gap studies ( Kinnunen et ai. 2004a). as well as density pro- 
file (Staiic e^^j2005a|^n^hermodvnarni^tudie^Kinast 
et ai, 2005) all appear to be compatible with this theory. 

The material in this Review is viewed as the first of many 
steps in a long process. It was felt, however, that some con- 
tinuity should be provided from another community which 
has addressed and helped to develop the BCS-BEC crossover, 
since the early 1990's when the early signs of the cuprate 
pseudogap were beginning to appear. As of this writing, it 
appears likely that the latest experiments on cold atoms probe 
a counterpart to this pseudogap regime. That is, on both sides, 
but near the resonance, fermions form in long lived metastable 
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pair states at higher temperatures than those at which they 
Bose condense. 
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APPENDIX A: Derivation of the T = variational 
conditions 



a) Z = 



b) 






FIG. 45 Self energy contributions (a) and response diagrams for 
the vertex correction corresponding to Epg (b). Heavy lines are for 
dressed, while light lines are for bare Green's functions. Wavy lines 



The ground state wavefunction is given by 

- e-^'/2+A6+ ^ jjkK + i;k4«tk)|0> (Al) 

Thus one has three variational parameters A, Wk, Wk with nor- 
malization condition + v^. = 1. We can define 6'k such 
that Uk — sin 6'k and Wk ~ cos ^k. so that we have only two 
variational parameters, A and ^k- 
We then have 



(*°|H-/i7V|*") = (iy-2Ai)A2 + 2^(ek-Ai)«; 



+ 2U^ UkWkWk'Wk' + 25A ^ -UkWk 

k,k' k 

where we have chosen i^k = 1 in Eq. il l\ for simplicity. 
Let us define 



k 

The variational conditions are 

5(H- fiN) 



SX 



and 



S{H- fiN) 
6K 



= 



(A2) 



(A3) 



(A4) 



These conditions yield two equations: a relationship between 
the two condensates 



ffA 1 g2 



Age U v — 2n 
and the T = gap equation: 

,2 



(A5) 



Similarly, we can write an equation for the total number of 
fermions 

= (i'°|2^6[6k + 2^aj,ak|«'°)=2A2 



E 



1 - 



Ck - M 



where 



(A7) 



with Asc = Asc — 5A. Thus we have obtained the three equa- 
tions which characterize the system at T = 0. It should be 
clear that the variational parameter A is the same as the con- 
densate weight so that A = (j)m = (&q=o)- 



APPENDIX B: Proof of Ward Identity Below Tc 

We now want to establish that the generalized Ward identity 
in Eq. J106> is correct for the superconducting state as well. 
The vertex 6Apg may be decomposed into Maki-Thompson 
(AIT) and two types of Aslamazov-Larkin {ALi, AL2) dia- 
grams, whose contribution to the response is shown here in 
Fig.Elh- We write 



SAr, 



5Kmt+5\\l + 5A\l{^). (B1) 



Using conventional diagrammatic rules one can see that the 
MT term has the same sign reversal as the anomalous super- 
conducting self energy diagram. This provides insight into 
Eq. ( I107> . Here, however, the pairs in question are noncon- 
densed and their internal dynamics (via tpg as distinguished 
from tsc) requires additional ALi and AL2 terms as well, 
which will ultimately be responsible for the absence of a 
Meissner contribution from this normal state self energy ef- 
fect. 
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Note that the AL2 diagram is specific to our GGa scheme, 
in which the field couples to the full G appearing in the T- 
matrix through a vertex A. It is important to distinguish the 
vertex A which appears in the AL2 diagram from the full EM 
vertex A-^*^. 

In particular, we have 



A + 5Ar 



(B2) 



where the sign change of the superconducting term (relative 
to A^^^) is a direct reflection of the sign in Eq. ( I107l i. 

We now show that there is a precise cancellation between 
the MT and ALi pseudogap diagrams at Q = 0. This can- 
cellation follows directly from the normal state Ward identity 



where are the Fourier components of the order parame- 
ter 5'(r,t), e = "^^J" , ^1 is the GL coherence length, fep ~ 
[l/T^7r^]|C(3), T* is the critical temperature when feedback 
effects from the quartic term are neglected, and N{0) is the 
density of states at the Fermi level in the normal state. The 
quantity (l^'qp) is determined self-consistently via 



(C2) 



where F\'^] is evalu^te^ji^^he_Jtotree_a£gro?^^ 

jic et al, 20 03bj . Our self consistency condition can be then 
written as 



Q-\{K,K + Q)^Go\K)-GqHK + Q), (B3) 
which implies 

Q ■ [5K\l{K, K + Q) + 6Amt{K, K + Q)] = Q (B4) 

so that 5A\j^{K^K) = —5Amt{K,K) is obtained exactly 
from the Q limit. 

Similarly, it can be shown that 

Q-K{K,K + Q) = G-^{K)-G-\K + Q) (B5) 

The above result can be used to infer a relation analogous to 
Eq. (IB4> for the AL2 diagram, leading to 



6Xpg{K,K) = -5AMT{K,K) 



(B6) 



which expresses this pseudogap contribution to the vertex 
entirely in terms of the Maki-Thompson diagram shown in 
Fig. 145b . It is evident that (5Aj\/t is simply the pseudogap 
counterpart of 5Asc, satisfying 



-5Kmt{K, K) = 



dk 



Therefore, one observes that for T < Tc 



5ApgiK,K) = 



dk 



(B7) 



(B8) 



which establishes that Eq. il06\ applies to the superconduct- 
ing phase as well. As expected, there is no direct Meissner 
contribution associated with the pseudogap self-energy. 



APPENDIX C: Quantitative Relation between BCS-BEC 
crossover and Hartree-TDGL 

In this appendix we make more precise the relation between 
Hartree-approximated TDGL theory and the T-matrix of our 
GGq theory. The Ginzburg-Landau (GL) free energy func- 
tional in momentum space is given by (Hassing and WiUcinsi 
Il973h 



yE*qi*q2*<l3*qi+q.^q3, (CI) 



T 

iV(0 



e + 6oE(l*q'ol')+^iV 
q' 



(C3) 



If we sum Eq. ( IC3> over q and identify X^qd^qop) with A^, 
we obtain a self-consistency equation for the energy "gap" (or 
pseudogap) A above Tc 



E 



T 



N{0) 



2„2 



e + boA^ + £,iq 



^(0) -A*pa.r(T)+e?g2 ' 



where 



(C4) 



(C5) 



(C6) 



Note that the critical temperature is renormalized downward 
with respect to T* and satisfies 



(C7) 



To compare with GL theory we expand our T-matrix equa- 
tions to first order in the self energy correction. The T-matrix 
can be written in terms of the attractive coupling constant U 

as 



t{Q) = 



U 



1 + UxoiQ) + USxiQ) ' 



(C8) 



where 



Defining 



Xo{Q) = J2^oiQ-K)Go{K) . (C9) 

K 



A2 = -^t(Q) 

Q 



(CIO) 



we arrive at the same equation as was derived in Section llLCl 
« -Go(-X)A2 (Cll) 
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where one can derive a self consistency condition on in 
terms of the quantity Sx{0) (which is first order in E), which 
satisfies 

5x{0) = -&oiV(0)A2, (C12) 

implying that 

which coincides with the condition derived earlier in Eq. (jC^)- 

APPENDIX D: Convention and Notation 
1. Notation 

We follow standard notations as much as possible. They are 
summarized below. 
Ep — Fermi energy 
kp — Fermi momentum 
h — Planck constant 
kB — Boltzmann constant 
c — Speed of light 
e — Electron charge 

Tc — Critical temperature for (superfluid/superconducting) 
phase transition 

T* — Pair formation or pseudogap onset temperature. 
T — Temperature 

fi, fipair, f^B — Fermionic, pair and bosonic chemical po- 
tential, respectively. 

A — Fermionic excitation gap 

Asc — Superconducting/superfluid order parameter 

Apg — Pseudogap 

Four vector K = (iw„, k), J^k = ^ J2n I^k' where uj^ = 
{2n + l)7rT is the odd (fermionic) Matsubara frequency. 

Four vector Q = (ifi^, q), J^q = T J2n I^q' where fi„ = 
2n7rT is the even (bosonic) Matsubara frequency. 

f{x) = l/(e^/'^^"^ + 1) — Fermi distribution function 
b{x) — l/{e^/^'^'^ — 1) — Bose distribution function 
G{K), Go{K) — Full and bare Green's functions for 
fermions 

D{Q), Dq{Q) — Full and bare Green's functions for Fes- 
hbach bosons. 

E(i4r), Y,b{Q) — Self energy of fermions and bosons, re- 
spectively. 

x{Q) — P^ir susceptibility 
t{Q) — T-matrix 

A, A — Bare and full vertex, respectively. 
I/, vq — Magnetic detuning parameter. 
Us, a* — S-wave inter-fermion scattering length 
U, Ua — Renormalized and bare strength of the attrac- 
tive interaction between fermions in the open-channel. Ua oc 

tts^V ~ +Oo). 

Uc — Critical coupling strength at which the two-body 
scattering length diverges. 

A = {(f), A) — Four vector potential 
J = (p, J) — Four current 



K^"{Q) — Electromagnetic response kernel 
ek = k'^ /2m — Bare fermion dispersion in free space 
= Ek — — Bare fermion dispersion measured from 

chemical potential. 

i?k — Bogoliubov quasiparticle dispersion 

= 5(1 + a/£^k), Ck/Ei,) — Coherence 

factor as given in BCS theory. 

71*°* = n + 2(nfc + rtf,o) — Total, fermion, finite momentum 

boson, condensed boson density, respectively. 
V{r) = imo^^r^ — Harmonic trap potential 
We use alphabetic ordering when multiple authors appear 

together 

We always refer to the absolute value when we refer to the 
interaction parameters U, Uq, g, go as increasing or decreas- 
ing. 

2. Convention for units 

Throughout this Review, we use the convention for units 
where it is not explicitly spelled out: 

h = kB = c = 1. 

Ep = Tp = kp = 2m = 1. 

In a harmonic trap, Ep, Tp and kp are defined by the non- 
interacting Fermi gas with the same total particle number N . 
Ep = fiaj{3Ny^^ . We further take the Thomas-Fermi radius 
Rtf — 1- 

In this convention, the fermion density in a homogeneous 
gas in three dimensions is n = I/Stt^, the harmonic trap 
frequency w = 2, and the total fermion number in a trap is 
N = 1/24. 

In the two-channel model, the units are Ep/kp for U and 
Uq, and Ep / k^^ for g and go. 

Our fermionic chemical potential /i is measured with re- 
spect to the bottom of the energy band, which leads to (i) 
ji = Ep in the non-interacting limit at T = 0, and (ii) /i 
changes sign when the system crosses the boundary between 
fermionic and bosonic regimes. 

3. Abbreviations 

BCS — Bardeen-Cooper-Schrieffer (theory) 

BEC — Bose-Einstein condensation 

GL — Ginzburg-Landau (theory) 

TDGL — Time-dependent Ginzburg-Landau (theory) 

RPA — Random phase approximation 

LDA — Local density approximation 

GP — Gross-Pitaevskii 

TF — Thomas-Fermi (approximation) 

NSR — Nozieres and Schmitt-Rink 

FLEX — FLuctuation EXchange 

AB — Anderson-Bogoliubov (mode) 

PG — Pseudogap 

SC — Superconductor 

AFM — Antiferromagnet 

FB — Feshbach boson 

LSCO — Lai_:rSr^Cu04 



46 



BSCCO — BiaSraCaCuzOs+a 
YBCO — YBasCuaOr-a 

ARPES — Angle-resolved photoemission spectroscopy 
STM — Scan tunneling microscopy 
RF - Radio frequency (spectroscopy) 
SIN — Superconductor-insulator-normal metal (tunneling 
junction) 

SI — Superconductor-insulator (transition) 
3D — Three dimensions 
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